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ABSTRACT 

This  report  summarizes  the  technical  effort  accomplished  in 
FY83  on  a  program  to  make  range-resolved  measurements  of  chemical 
agents  or  pollutants  distributed  within  the  atmosphere.  A 
description  of  the  modified  receiver  electronics  which  permit 
rapid  accumulation  of  a  very  large  number  of  signal  returns  is 
described.  The  performance  of  a  rapidly  tunable  CO^7 laser  is 

l 

discussed.  Theoretical  analyses  of  coherent  laser  radar  systems 
including  mixing  efficiencies,  speckle  effects,  and  waveform 
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INTRODUCTION 


This  is  the  FY83  final  report  on  a  program  entitled 
"Coherent  Laser  Radar  Remote  Sensing"  supported  by  the  Air  Force 
Engineering  and  Services  Center  ( AFESC)  located  at  Tyndall  Air 
Force  Base.  This  effort  is  part  of  a  program  to  assess  the 
utility  of  using  a  coherent  laser  radar  for  the  range-resolved 
detection  and  identification  of  chemical  agents  in  the 
atmosphere.  This  work  was  begun  in  April  1982  and  will  continue 
through  FY84.  Previous  research  has  been  documented  in  the  FY82 
final  report1 . 

Specific  tasks  required  under  this  program  include 
1)  modifications  to  the  existing  transportable  laser  radar 
receiving  electronics  to  permit  recording  of  aerosol  backscatter 
in  several  range  bins  on  a  single  pulse,  2)  field  measurement  of 
atmospheric  aerosol  backscatter  cross-sections  at  single 
wavelength  near  10.6um  with  the  modified  electronics  package, 

3)  modification  of  the  laser  transmitter  to  permit  rapid 
switching  between  at  least  two  separate  wavelengths. 

In  section  II  of  this  report  we  describe  the  changes  to  the 
receiver  electronics  that  are  required  and  report  on  their 
status.  Section  III  describes  the  progress  made  in  constructing 
a  laser  rapidly  tunable  in  wavelength  by  using  an  intracavity 
etalon.  Several  theoretical  issues  have  been  investigated  and 
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are  summarized  in  Section  IV.  These  issues  include  optimization 
of  heterodyne  efficiency,  decorrelation  effects,  and  the  use  of 
other  than  a  pulse  waveform  for  achieving  range  resolution. 
Detailed  technical  descriptions  of  these  issues  have  been 
included  in  three  appendices. 

II .  MODIFICATIONS  TO  RECEIVER  ELECTRONICS 

The  rationale  for  using  a  compact,  coherent  CC2  laser  radar 
to  remotely  sense  concentrations  of  chemical  agents  was 
established  in  a  previous  report1 .  The  basic  radar  performance 
theory  for  range-resolved  aerosol  return  measurements  was 
presented  therein,  along  with  experimental  verifications  of  the 
key  premises  of  that  theory:  namely,  signal  strengths  and 
statistics  agreed  well  with  theory.  However,  rapid  data 
collection  was  hampered  by  the  fact  that  the  original 
electronics  had  been  designed  to  operate  on  hard  target 
returns.  When  operating  against  solid  targets  it  is  sufficient 
to  record  the  strength  and  time  of  return  of  the  single  largest 
signal  within  the  range-gate  window.  For  the  aerosol  problem, 
the  optical  radar  should  be  operated  much  like  a  conventional 
microwave  radar,  in  which  the  signal  propagates  through  space 
encountering  many  targets.  Thus  samples  of  the  return  signal 
should  be  recorded  at  many  times  from  a  single  pulse. 

The  initial  data  were  recorded  in  this  mode  through  the  use 
of  a  commercial  digital  oscilloscope  with  floppy-disk  memory. 
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This  set-up  permitted  recording  signals  at  75  different  ranges 
between  0.2km  and  1.5km  on  each  pulse.  Although  the  system  runs 
at  a  pulse  repetition  frequency  ( PRF )  of  20kHz,  the  collection 
time  for  2,500  samples  (pulses)  of  data  was  approximately  3 
minutes  in  this  mode  owing  to  both  a  memory  dump  rate  limitation 
and  a  memory  size  limitation  in  the  oscilloscope.  Receiver 
electronics  have  been  built  which  will  increase  both  the  data 
acquisition  rate  and  the  amount  of  data  stored.  The  system  will 
run  at  an  8.66kHz  pulse  rate  and  collect  and  record  signal 
strengths  from  60  samples  spaced  at  100ft  (30.5m)  range 
intervals  on  each  pulse.  The  acquisition  window  of  6000ft 
(1800m)  may  be  placed  anywhere  within'  a  setting  of  0  and 
20,000ft  (6,1km).  In  practice  returns  from  ranges  less  than 
300m  are  ignored  by  the  system. 

This  set  of  electronics  is  complete  and  is  awaiting  final 
"debugging"  checks.  When  operational  this  system  will  permit 
acquiring  data  in  times  of  approximately  1  second,  which,  when 
averaged,  should  permit  making  measurements  with  an  accuracy  of 
1%.  The  experiments  on  the  aerosol  backscatter  cross-section 
will  then  be  repeated,  and  statistics  gathered  on  the 
variability  of  this  cross-section  at  a  given  time  but  at 
different  look  angles  and  also  as  a  function  of  different 
weather  conditions. 

Single  v.’ave  length  absorption  measurements  will  then  be 
made.  A  large  absorption  cell  will  be  placed  at  a  range  of 
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approximately  3/4km  from  the  radar.  With  the  cell  empty  refer¬ 
ence  measurements  of  the  backscatter  will  be  made.  By  filling 
the  cell  with  an  adsorbing  gas  (such  as  ethylene  or  sulfur  hexa¬ 
fluoride),  the  sensitivity  of  the  system  may  be  established. 

This  technique  simulates  dual  wavelength  Differential  Absorption 
Lidar  (DIAL)  experiments.  If  the  issues  addressed  in  section  ill 
of  this  report  can  be  solved  this  year  a  full  dual  wavelength, 
range  resolved,  remote  sensing  experiment  will  be  attempted, 
in.  FREQUENCY  AGILE  CO?  LASER 

During  FY83,  we  have  implemented  the  intracavity  tunable 
etalon  technique  for  rapid  tuning  of  the  C02  laser  between 
transitions.  We  employed  a  commercial  etalon  (Burleigh  TL15IR) , 
with  a  laboratory  testbed  laser.  The  etalon  mirrors  are  coat. id 
for  80%  reflectance  at  10. 6pm  and  spaced  for  a  free  spectral 
range  of  5C  cm"1  .  This  allows  tuning  over  the  entire  p  branch  of 
the  10um  band. 

The  optical  configuration  is  shown  schematically  in  Figure 
1.  As  the  etalon  is  *-uned  monotonically,  the  laser  output  shifts 
from  line  to  line,  mostly  between  successive  lines  of  the  10um  p 
branch.  In  the  low-gain  region  of  this  spectrum,  however  the 
laser  will  occasionally  emit  at  R-branch  wavelengths  or  even  in 
the  9ym  band.  This  behavior  is  repeatable,  as  was  demonstrated 
by  applying  a  ramp  voltage  to  the  etalon.  The  non-p  branch 
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oscillation  could  be  suppressed,  if  desired,  by  appropriate 
coating  of  the  laser  output  mirror,  ensuring  that  all  such 
competing  lines  are  below  oscillation  threshold.  Alternatively, 
the  etalon  spacing  could  be  decreased  by  approximately  a  factor 
of  4  to  bring  these  lines  within  the  free  spectral  range.  This, 
however,  would  also  reduce  the  spectral  resolution 
proportionally  and  hence  the  reliability  with  which  a  given  line 
could  be  accessed. 

If  this  type  of  laser  is  to  be  used  in  a  differential 
absorption  LIDAR  (DIAL)  system,  it  will  have  to  be  capable  of 
switching  rapidly  between  at  least  two  prescribed  transitions 
(on-  and  off-  resonance).  To  demonstrate  the  capability,  the 
laser  was  tuned  to  P(2G),  and  a  squarewave  input  pulse  of 
variable  amplitude  was  applied  to  the  DC  amplifier  used  to  drive 
the  etalon.  By  this  means,  we  were  able  to  switch  back  and 
forth  between  P(20)  and  any  of  17  other  lines  at  rapid  rates, 
limited  presumably  only  by  the  response  of  the  piezoelectric 
transducer.  Included  in  this  group  was  the  pair  P(20),  P(14) 
which  can  be  used  for  DIAL  measurements  of  C2H4  (ethylene) 
concentration.  Another  possible  pair  for  C2H4  detection  is 
P(12),  P(14)  for  which  rapid  line  shifting  was  also  observed. 

The  primary  advantage  of  this  electronic  tuning  technique 
over  conventional  mechanical  methods  is  its  speed.  To  quantify 
this,  a  Ge:Au  detector  was  used  to  monitor  the  transient  output 
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as  the  laser  was  tuned  from  P(14)  to  P(20).  As  shown  in  Figure 
2 ,  the  amplitude  fluctuations  decayed  in  several  milliseconds. 

No  attempt  was  made  to  damp  these  fluctuations. 

One  problem  encountered  with  this  approach  is  that  the 
cavity  length  is,  in  general,  not  resonant  for  both 
transitions.  Accordingly,  the  maximum  power  cannot  be  obtained 
from  each  of  the  sequentially  tuned  lines  without  also  adjusting 
the  total  cavity  length.  We  addressed  this  problem  by  applying 
the  squarewave  input  to  the  D.C.  amplifier  controlling  the  main 
cavity  length  (as  well  as  to  the  etalon  driver)  and  adjusting 
the  gain  to  maximize  the  laser  power. 

Another  approach,  which  we  have  also  implemented,  is  to 
change  the  gross  cavity  lengtn  such  that  the  two  relevant 
transitions  may  be  simultaneously  resonant.  For  the  wavelengths 
involved,  the  required  length  change  is  invariably  less  than  a 
centimeter,  so  that  we  were  able  to  accomplish  it  by  mounting  a 
cavity  mirror  on  a  piezoeletric  inchworm  translator  (1-inch 
total  travel).  This  cavity  retro-fit  enabled  us  to  tune 
sequentially  between  each  member  of  any  predetermined  accessible 
line  pair  without  also  adjusting  the  cavity  length  with  each 
line  change. 

Another  concern  when  intracavity  etalons  are  used  is 
thermal  stability.  Indeed  the  laser  required  approximately  a 
half-hour  warm-up  before  approaching  thermal  equilibrium.  After 


this  time,  the  etalon  required  periodic  bias  adjustment  to 
maintain  operation  on  a  given  line  pair,  but  the  required  step 
voltage  remained  constant  and  was  repeatable  on  a  day-to-day 
basis.  It  is  thus  possible,  once  operation  on  a  given  line  is 
achieved,  to  switch  predictably  to  any  other  accessible 
transition.  The  required  bias  control  would  be  achieved  in  an 
operational  system  by  a  servo  control  loop. 

One  disadvantage  of  the  etalon  technique  of  wavelength 
tuning  is  the  large  insertion  loss  of  the  etalon,  which  we 
measured  to  be  25%  per  pass.  This  resulted  in  a  decrease  in 
laser  power  from  approximately  7  watts  to  2  watts.  That  this 
power  loss  is  so  large  is  due  to  the  fact  that  the  testbed  laser 
is  relatively  short  (active  length  «  50cm),  and  with  the 
intracavity  etalon,  it  operates  very  near  threshold.  We  are 
currently  fabricating  a  longer  laser  (active  length  «  120cm)  to 
be  used  as  a  testbed,  for  which  this  power  reduction  should  be 
much  less. 

Another  approach  for  reducing  the  etalon  insertion  loss  is 
to  reduce  the  mirror  reflectance.  The  resulting  decrease  in 
power  dissipation  by  the  etalon  should  also  reduce  any  thermal 
distortion.  This  change  would  also  reduce  the  etalon's  finesse 
and  hence  its  resolving  power,  but  a  modest  reduction  in 
reflectance  might  improve  the  overall  performance. 


in  conclusion,  we  have  demonstrated  the  intracavity  etalon 
technique  for  rapid  wavelength  tuning  of  a  CO2  laser.  Areas  for 
further  study  include  multiple  wavelength  operation  using 
multiple  step  driver  waveforms  and  the  effect  of  etalon  mirror 
reflectance  on  laser  performance.  In  addition,  there  remains 
the  formidable  task  of  frequency  locking  a  local  oscillator  to 
the  frequency-agile  transmitter.  T.n  this  regard,  the  pulsed-cw 
waveform  for  pulsed  homodyne  operation  might  be  the  most 
practical  solution, 
iv.  theoretical  analyses 

The  theoretical  effort  during  FY83  was  directed  toward 
design-optimization  issues  for  aerosol  measurement  laser 
radars.  In  particular,  because  of  the  weakness  of  aerosol 
returns  relative  to  hard-target  returns,  against  which  compact 
CO2  laser  radars  are  customarily  employed,  maximization  of 
heterodyne-mixing  efficiency  is  critical  for  chemical  agenc 
detection.  Moreover,  because  fielded  systems  may  have  to  be 
multi-functional,  mixing  efficiency  should  be  considered  in 
conjunction  with  the  spatial  resolution  of  the  radar,  a 
characteristic  critical  to  its  hard-target  imager  function. 
Furthermore,  Doppler-imaging  radars,  and  range-imaging  radars 
which  use  high  time-bandwidth  (TW)  product  waveforms,  may  also 
be  part  of  the  multi-functional  panoply.  Thus,  the  effects  of 


aerosol  decorrelation  on  these  long  dwell  time  systems  should  be 
understood,  so  that  the  potential  of  such  radars  as  chemical 
detectors  can  be  assessed. 

Results  have  been  obtained  and  presented  in  all  of  the 
theoretical  areas  described  above3-5.  Their  salient  features 
are  described  below;  the  full  papers  cited  appear  as  appendices 
to  this  report. 

A.  Antenna  Patterns  and  Mixing  Efficiencies3 

The  choice  of  transmitter  and  local-oscillator  beam 
patterns  affects  both  the  spatial  resolution  capability  (through 
the  transmitting  and  receiving  antenna  patterns)  and  the 
carr ier-to-noise  ratio  (through  the  heterodyne-mixing 
efficiency)  of  a  coherent  laser  radar.  We  have  made  the  first 
coordinated  assessment  of  these  effects.  A  monostatic 
shared-optics  radar  employing  an  unobscured  circular-pupil 
objective  (diameter  D),  a  circular-pupil  detector,  and  Gaussian 
transmitter  and  local-oscillator  beams  was  considered.  In 
far-field  operation,  parameter  values  were  obtained  yielding 
mixing  efficiencies  in  excess  of  80%  and  25%  for  glint  targets 
and  speckle  targets,  respectively,  and  better  than  0.8  X/D 
angular  resolution,  where  X  is  the  laser  wavelength. 

B .  Resolution  and  Carr ier-to-Noise  Ratio4 , 5 

The  long  pixel  dwell  times  of  Doppler  laser  radars  and  high 
TW  product  range-imaging  laser  radars  make  them  susceptible  to 


9 


performance  degradation  arising  from  target  decorrelation.  We 
have  included  target  decorrelation,  in  addition  to  speckle 
effects,  in  analyses  of  velocity  and  range  resolution  and 
carrier-to-noise  ratio  (CNR)  for  such  radars.  Of  particular 
interest  for  chemical  agent  detection  is  the  fact  that 
chirped-pulse  range  resolution  degrades  rapidly  when  the  chirp 
duration  T  exceeds  the  aerosol  correlation  time  (typically  about 
lus),  but  the  onset  of  significant  CNR  loss  does  not  occur  until 
this  correlation  time  becomes  smaller  than  the  reciprocal  chirp 
bandwidth  1/W.  Because  high  TW  operation  with  T  >  lus  and  W  > 
10MHz  is  assumed,  chirped-pulse  aerosol  measurements  are  likely 
to  suffer  resolution  loss  out  not  CNR  loss. 

V.  SUMMARY  AND  RECOMMENDATIONS 

In  FY82  range-resolved  signal  returns  from  naturally 
occurring  atmospheric  aerosols  were  observed  from  ranges  up  to 
1.5km  with  the  transportable  coherent  CO2  laser  radar  system. 

The  magnitude  of  the  signals  was  approximately  that  expected 
theoretically.  In  FY83  modifications  to  the  receiver 
electronics  were  designed  and  constructed  that  will  permit 
accumulating  range-resolved  data  at  an  8.7kHz  rate.  In  addition 
an  intracavity  etalon  technique  was  used  to  demonstrate  rapid 
tuning  between  pairs  of  lines  across  the  9ura  to  Hum  spectral 
band . 

Theoretical  analyses  of  waveforms  other  than  pulse  show 
decorrelation  effects  caused  by  the  motion  of  the  aerosols 
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during  the  measurement  interval.  These  decorrelations  produce  a 
loss  of  range  resolution  but  no  loss  in  carr ier-to-noise  ratio. 

The  optimum  ratio  of  the  relative  sizes  of  the  transmitter  and 
receiver  apertures  was  also  obtained  by  analysis. 

Work  in  FY84  is  aimed  at  assessing  the  variability  of  the 
aerosol  cross-sections  and  demonstrating  a  differential  absorption 
measurement  at  a  single  wavelength  of  a  known  gas  released  into 
the  atmosphere.  If  the  problem  of  having  the  local  oscillator 
wavelength  track  the  transmitter  wavelength  can  be  solved  in  a 
timely  fashion,  then  dual  wavelength  absorption  measurements  will 
be  made. 
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ABSTRACT 

Transmitting  and  receiving  antenna  pattern  expressions 
are  developed  for  a  monostatic  coherent  laser  radar.  These 
patterns  are  rigorously  related  to  the  standard  glint-target 
and  speckle-target  monostatic  radar  equations  to  yield 
formulas  for  glint  and  speckle  target  mixing  efficiencies. 
Numerical  calculations  for  spatial  resolution  and  mixing 
efficiency  are  presented  for  far-field  operation  using  an 
untruncated-Gaussian  transmitter  beam  and  a  Gaussian  local- 
oscillator  beam  that  has  suffered  both  photodetector  and 
objective-pupil  truncation.  Closed-form  expressions  for 
spatial  resolution  and  mixing  efficiency  for  an  all  Gaussian 
baam  system  are  employed  io  assess  the  need  for  dynamic 
focusing  in  a  large-aperture  radar  operating  in  the  near 
field. 


1 . 0  INTRODUCTION 

Future  requirements  for  high-resolution  multifunction  coherent  C02  laser 
radars  dictate  a  need  for  accurate  system  design  methodologies.  In  particular, 
the  choice  of  transmitter  and  local  oscillator  spatial  beam  patterns  directly 
affects  both  the  radar's  spatial  resolution  capability  (through  the  transmitting 
and  receiving  antenna  patterns)  and  its  carr ier-to-noise  ratio  (through  the 
heterodyne-mixing  efficiency).  Previous  studies  have  not  made  a  coordinated 
assessment  of  these  effects.  There  is  published  work1, z, 3  on  mixing 


1.  Degnan,  J.J.,  and  Klein,  B.J.,  "Optical  Antenna  Gain  2:  Receiving  Antennas," 
Appl.  Opt.  13,  2397-2401  (1974). 

2.  Fink,  D. ,  "Coherent  Detection  Signal-to-Noise,"  Appl.  Opt.  14,  689-690 

(1975).  — 

3.  saga,  N.,  Tanaka,  K.,  and  e'ukumitsu,  u. ,  "Diffraction  ot  a  Gaussian  Beam 
through  a  Finite  Aperture  Lens  and  the  Resulting  Heterodyne  Efficiency," 

Appl'.  Opt.  20,  2827-2831  (  1981). 

"This  work  was  sponsored  bv  the  Department  of  the  Air  Force. 

"The  US  Government  assumes  no  responsibility  for  the  information  oresented". 


A-1 


efficiency  as  a  function  of  local-osci llator  beam  pattern  and  photodetector 
diameter,  but  this  work  neglects  the  effects  cf  target  speckle  and  transmitter 
beam  pattern  on  mixing  efficiency,  and  those  of  local-oscillator  beam  pattern  and 
photodetector  diameter  on  the  receiving  antenna  pattern.  The  work  that  has  been 
published1*,5,0  on  mixing  efficiency  including  transmitter  and  speckle  effects 
often  assumes  untruncated  Gaussian  beams,  and  generally  does  not  examine  the 
radar's  spatial  resolution.  This  paper  addresses  the  joint  dependence  of  spatial 
resolution  and  mixing  efficiency  on  the  transmitter  and  local  oscillator  beam 
patterns,  eliminating  some  of  the  gaps  left  by  earlier  treatments. 

Section  2.0  presents  the  system  structure  to  be  examined  and  the  fundamental 
assumptions  that  are  made.  The  analysis  begins,  in  Section  3.0,  with  the  devel¬ 
opment  of  antenna-pattern  expressions.  These  are  then  rigorously  related  to  the 
standard  glint-target  and  speckle-target  monostatic  radar  equations  to  yield  for¬ 
mulas  for  glint  and  speckle  target  mixing  efficiencies.  Numerical  calculations 
are  presented  in  Section  4.0  for  spatial  resolution  and  mixing  efficiency  for 
far  field  operation  using  an  untruncated-Gaussian  transmitter  beam  and  a  Gaussian 
local-oscillator  beam  that  has  suffered  both  photodetector  and  objective-pupil 
truncation.  Closed-form  expressions  for  spatial  resolution  and  mixing  efficiency 
for  an  all  Gaussian  beam  system  are  employed  in  Section  5.0  to  assess  the  need 
for  dynamic  focusing  in  a  large-aperture  radar  operating  in  the  near  field. 

2.0  SYSTEM  STRUCTURE 

The  system  structure  we  shall  consider,  shown  schematically  in  Fig.  1,  is  a 
monostatic  shared-optics  coherent  laser  radar  using  a  single  photodetector  and 


4.  An  all  Gaussian  oeam  signa 1-to-noi se  ratio  calculation  for  a  nonostatic  radar 
viewing  a  speckle  target  appears  in  Sonnensche in ,  C.M.,  and  Horrigan,  F.A., 

"S ignal-to-Noise  Relationships  for  Coaxial  Systems  that  Heterodyne  back- 
scatter  from  the  Atmosphere,"  Appl.  Opt.  1_0,  1600-1604  (  1971  ).  This  work 
neglects  atmospheric  turbulence,  and  does  not  examine  spatial  resolution. 

5.  All  Gaussian  beam  speckle-target  mixing  efficiency  including  the  effects  of 
atmospheric  turbulenge  was  treated  by  Yura,  H.T.,  "Signal-to-Noise  Ratio  of  . 
Heterodyne  Lidar  Systems  in  the  Presence  of  Atmospheric  Turbulence,"  Optica 
Acta  2_6,  627-644  (1979),  and  Clifford,  S.F.,  and  Wandzura,  S.,  "Monostatic 
Heterodyne  Lidar  Performance:  The  Effect  of  the  Turbulent  Atmosphere," 

Appl.  Opt.  2_0,  514-516  (1981).  Neither  work  studied  spatial  resolution. 

6.  Mixing  efficiency  for  truncated  Gaussian  beam  radars  has  been  analyzed  in 
Rye,  B.J.,  "Antenna  Patterns  for  Incoherent  Backscatter  Heterodyne  Lidar," 
Appl.  Opt.  18^,  1390-1398  (  1979)  and  calculations  have  been  reported  by 
Leader,  J.C.,  "Detection  Efficiency  for  Large-Aperture  Coherent  Laser 
Radars,"  Proc.  SPIE  300 ,  74-85  (1981  ),  and  Wang,  J.Y.,  "Heterodyne  Laser 
Radar  SNR  from  a  Diffuse  Target  Containing  Multiple  Glints,"  Appl.  Opt.  2 1 , 
464-475  (1982).  The  latter  two  studies  include  atmospheric  turbulence;  none 
of  these  papers  coraiders  spatial  resolution  per  se. 


The  transmitter  and  local-oscil- 


unobscured  circular  pupil  refractive  optics7, 
lator  beams  entering  the  transmit/receive  (TR)  switch  are  both  TEMqu  untrun¬ 
cated  Caussians.  These  beams  are  such  that  the  normalized  transmitter  spatial 
beam  pattern  as  it  leaves  the  objective  lens  is 

£t(p>  *  (2/»aT2)1/2  exp[-| o  j  2  (a^2  +  jk/2R)],  (li 

and  the  normalized  local-oscillator  spatial  beam  pattern  as  it  impinges  on  the 
photodetector  is 

£l<p)  -  (2/t3dl2)1/2  exp(-j  e  j  2  (a^2  -  jk/2i  )  ] ,  (2) 

for  e  ■  (x,y)  the  transverse  coordinate  vector,  and  k  *  2»/X  the  wavenumber  at 
the  radar  wavelength  X.  Here,  a-p  and  a^  are  e~2  intensity  radii  of  the 
transmitter  and  local  oscillator  beams  at  the  objective  lens  and  the  photodetec¬ 
tor,  respectively.  The  transmitter  beam  (1)  ha3  a  phase  radius  of  curvature  cor¬ 
responding  to  a  geometric-optics  focus  located  R  (R>0)  meters  in  front  of  the 
objective,  implying  a  Gaussian  beam  waist  located  R/(l+Q~2)  meters  in  front  of 
the  objective  where  a  «  ka-j2/2R.  The  local -oscillator  beam  (2)  has  a  phase 
radius  of  curvature  such  that  were  backpropagated  through  the  optical 

system  the  beam  exitinq  the  objective  lens  would  have  e-2  intensity  radius 

aL  *  *4/»aDL'  (3) 

and  the  same  phase  radius  of  curvature 

R  -  (f_1  -  t"1)"1,  (4) 

as  the  transmitter  beam,  where  f  is  the  focal  length  of  the  objective  lens. 

The  preceding  Gaussian-beam  descriptions  ignore  the  truncation  effects 
imposed  by  the  finite  photodetector  diameter  d  (on  the  local  oscillator),  and  by 
the  finite  pupil  diameter  D  (on  the  transmitter  and  the  backpropagated  local 
oscillator).  Because  transmitter  truncation  entails  a  loss  of  transmitter  power, 
which  is  generally  at  a  premium,  we  shall  assume  a  condition  that 

ensures  transmitter  truncation  is  negligible  in  terms  of  both  power  and  target- 
plane  beam  pattern.  There  is  no  similarly  high  premium  on  the  local-oscillator 
power,  so  in  general  that  beam  may  be  subject  to  both  of  the  truncation  effects 
noted  above.  In  what  follows,  (3)  will  be  used  to  define  a and  (4)  to  relate 


7.  The  formulas  we  derive  apply  as  well  to  unobscured  reflective  optics  such  as 
those  employed  in  the  Lincoln  Laboratory  test  bed  radar,  see  Harney,  R.C., 
and  Hull,  R.J.,  "Compact  Infrared  Radar  Technology,"  Proc.  SPIE  227,  162-170 
(1980). 
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R ,  t,  and  the  objective  focal  length  f,  even  though  the  backpropagated 
local-oscillator  field  leaving  the  objective  pupil  need  not  be  Gaussian.  Thus, 
the  parameters  at,  and  R  cease  having  the  physical  interpretations  given  with 
(3),  (4)  unless  a  _>  4apb  and  D  2  4al  both  prevail. 

In  any  case,  the  target  of  interest  will  be  assumed  to  De  a  stationary 
object  that  is  unresolved  in  range  (planar)  and  at  a  distance  L  m  in  front  of  the 
radar.  The  target  will  either  be  a  spatially  unresolved  on-axis  glint  object  of 
radar  cross-section  a  m2,  or  a  spatially  resclved  speckle  ODject  of  uniform 
diffuse  reflectivity  p  sr-^.  The  atmosphere  will  be  assumed  to  be 
non-turbu lent  with  an  extinction  coefficient  a  m-^.  All  losses  in  the  radar 
optics  (not  associated  with  local-oscillator  truncation)  will  be  neglected.  The 
photodetector  will  be  taken  to  have  quantum  efficiency  n. 

3.0  BASIC  RESULTS 

The  normalized  intermediate-frequency  (IF)  photocurrent  for  the  radar  shown 
in  Fig.  1  has  complex  envelope8 

r(t)»y+n(t)  (5) 

consisting  of  a  target  return  y  plus  local-oscillator  shot  nosie  n(  t ) .  Via  the 
antenna  theorem  for  heterodyne  detection9,  the  former  is  conveniently  expressed 
as  a  target-plane  integral8,10,11 

y  -  P r1/2J  dp  T(o)  £T'(?)  S.r*  (o  )  ,  (6) 

in  terms  of  the  transmitter  power  Pj,  the  transmitter  pattern  ' ( p  ) 
resulting  from  propagating  £x(p)  over  the  L-meter  atmospheric  path,  the 
receiver  pattern  §_r'(p)  resulting  from  truncating  and  renormalizing  5j/(p)  on 
the  detector  pupil  then  backpropagating  through  the  optical  system  and  over  the 
L-meter  atmospheric  path,  and  the  target's  complex-field  reflection  coefficient 


8.  The  svbcera  model  we  shall  employ  parallels  that  developed  in  Shapiro,  J.H., 
Capron,  B.A.,  and  Harney,  R.C.,  "Imaging  and  Target  Detection  with  a  Hetero¬ 
dyne-Reception  Optical  Radar,"  Appl.  Opt.  _2U,  3292-3313  (  1981  ). 

9.  Siegman,  A.E.,  "The  Antenna  Properties  ot  optical  Heterodyne  Receivers," 
Proc.  IEEE  54,  1350-1356  (1966). 

10.  Rye,  B.J.,  "Ref ract i ve- Index  Turbulence  Contribution  to  Incoherent 
Backscatter  heterodyne  Lidar  Returns,"  J.  Opt.  Soc.  Am.  71 ,  687-691  (1981). 

11.  Papurt,  D.M.,  Shapiro,  J.H.,  and  Harney,  R.C.,  "Atmospheric  Propagation 
Effects  on  Coherent  Laser  Radars, “  Proc.  SPIE  300 ,  86-99  (1981  ). 
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T(p).  Via  well  known  results  for  heterodyne  detection  ,  the  noise  ri(t)  is  a 
zero-mean  circulo-complex  white  Gaussian  process  of  spectral  height  hvQ/n,  in 
terms  of  the  photon  energy  at  the  radar's  optical  frequency  \>Q  »  c/X. 

3.1  GLINT  TARGET  RESULTS 


The  glint  target  under  consideration  is  spatially  unresolved,  thus  the  only 
performance  measure  to  be  evaluated  for  this  case  is  the  carrier-to-noise  ratio 
(CNR),  namely  tne  ratio  of  target-return  power  to  average  noise  power  in  a  S  Hz 
IF  bandwidth.  By  means  of  previous  work13,  we  know  that  the  CNR  obeys  the 
monostatic  radar  equation 

hPt  gt  oAr9  2  ^ 


CNR 


hv  B  4iL  4*L 
o 


(7) 


where 

Gt  -  4*L2  eaL  j  ST ' ( 3 ) | 2  .  (9) 

is  the  transmitter  antenna  gain,  and 


AR9  -  (iL)2  e°L  Ur'<5)|2  (9) 

is  the  effective  heterodyne-detection  glint-target  mixing  area.  Equation  (7)  i3 
the  standard  glint-target  result  if  Ar9  is  set  equal  to  the  objective  pupil 
area  »D2/4.  It  is  a  simple  Fourier  optics  calculation  to  show  that  Ar9  £ 
i0‘/4,  so  that  we  may  usefully  define 

ehet9  “  4ar9/*d2'  (10» 

as  the  glint-target  monostatic  radar  heterodyne  mixing  efficiency. 

3.2  SPECKLE  TARGET  RESULTS 

For  a  resolved  speckle  target  of  non-uniform  reflectivity  tho  average  target 
return  power  obeys8 'H 

<|y|2>  -  PTX2/  <3p  T(p )  |  £T  •  (  p  )  |  2  |£r*  <p  >  |  2  »  (ID 

where  7"  ( p )  is  the  average  intensity  reflection  profile  of  the  target.  Evidently, 

the  spatial  resolution  of  the  radar  (at  range  L)  is  determined  by  the  transmitter 


12.  Gagliardi,  R.  M. ,  and  Karp,  S.,  Optical  Communications  (Wiley,  New  York, 
1976)  Chap.  6. 

13.  Shapiro,  J.H.,  "Target-Reflectivity  Theory  for  Coherent  Laser  Radars," 
Appl.  Opt.  21,  3398-3407  (  1982  ). 
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antenna  pattern 


Ct(p  >  *  | £T * ( d )/£t * ( 0 ) | 2  ,  (12) 

and  the  receiver  antenna  pattern 

;R(P>  =  |iR'(J)/iR’(0)|2  ,  (13) 

where  we  have  assumed  that  the  denominator  terms  in  (12)  and  (13)  are  non-zero. 

In  particular,  the  nominal  sp  tial  resolution  of  the  radar  will  be  taken  to  equal 
the  j  o |  value  necessary  to  drive  the  composite  antenna  pattern  cT(p)cR(p) 
from  its  on-axis  value  of  unity  to  a  value  of  e-2.  (Note  that  our  system 
structure  forces  the  antenna  patterns  to  depend  only  on  |  p  J .  )  In  far-field 
operation  the  spatial  resolution  will  be  proportional  to  L,  i.e.,  the  radar  has  a 
fixed  far-field  angular  resolution. 

Suppose  that  the  target  profile  T(p)  corresponds  to  a  uniform  diffuse 
reflectivity  p,  viz.13  T  =  p/n,  over  the  spatial  resolution  cell  probed  by  the 
radar.  It  then  follows  that  the  speckle-target  CNR  obeys 


where 


CNR  » 


nP, 


hv  B 
o 


-2aL 


(14) 


Ars  =  (XL)2  e2aL  /  dp  |  §.T  ’  (  p  )  j  2  |iR’(p) 


(15) 


is  the  effective  heterodyne-detection  speckle-target  mixing  area, 
this  is  the  standard  CNR  result  if  Aps  *  *D2/4,  and  Fourier  optics 
used  to  prove  ARS  _<  »D2/4.  Thus,  we  define 


Once  again, 
can  be 


ehet3  3  4arS/»d2  '  <16> 

as  the  speckle-target  monostatic  radar  heterodyne  mixing  efficiency. 

4.0  FAR-FIELD  CALCULATIONS 

In  this  section  we  shall  present  calculations  for  the  case  of  far-field 
operation  (kD2/4L  <<  1)  with  the  radar  focused  at  infinity  (R  *  • ,  l  *  f). 

Here  our  main  effort  will  be  to  study  the  effects  of  local  oscillator  truncation 
on  spatial  resolution  and  mixing  efficiency. 

4.1  GLINT  TARGET  CALCULATIONS 

Under  the  propagation  conditions  cited  above,  the  basic  glint-target  results 
take  the  form  given  in  Section  3.1  with 


G 


T 


(17) 
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and  (after  some  manipulation) 


where  Jj  is  a  Bessel  function, 


and 


WL  *  4aL/D  ' 

y  -  Dd/2X  f  . 


(IB) 


(19) 


(20) 


The  dimensionless  parameters  ul  and  T  are  k®y  physical  entities:  ul  measures 
the  backpropagated  local  oscillator  size  (in  the  absence  of  photodetector 
truncation,  4aoL/di1>  relative  to  the  objective  diameter;  y  measures  the 
detector  size  relative  to  the  focal-plane  diffraction  width  of  the  objective. 

Both  parameters  can  be  expected  to  be  of  order  of  magnitude  unity  in  realistic 
system  designs. 

Two  important  special  cases  of  Eq.  (18)  are  easily  developed.  When  4aoL/d 
>>  1  (equivalent  to  8/*uly  >;>  U  the  local-oscillator  beam  (2)  has  an 
essentially  uniform  intensity  over  the  diameter  d  detector;  it  can  then  be  shown 
that 

Lhet9  *  t2<1  "  J0(*t))/»t]2  .  (21) 

On  the  other  hand,  when  4aQ[/d  £  1  and  4a[/D  £  1  (equivalent  to  8/tuly  <,  1 
and  ul  £  1)  there  is  no  truncation  of  the  local  oscillator  on  the  detector  nor 
is  there  any  (after  backpropagation )  at  the  objective  lens;  for  this  case 


In  Fig.  2,  we  have  plotted  £hetJ  !from  Eq.  (20)!  vs.  y  for  a  variety  of 
values.  For  a  fixed  transmitter  beam  and  objective  lens,  this  amounts  to 
plotting  glint-target  mixing  efficiency  vs.  normalized  photodetector  diameter  for 
various  normalized  local-oscillator  intensity  radii.  We  see  from  Fig.  2  that  for 
each  ul  value  there  is  an  optimum  y  value  at  which  c^et^  is  maximized,  we 
also  see  ehet9  1  0,8  is  attainable  within  the  parameter  set  we  have  plotted. 
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Note  that  as  the  photodetector  diameter  tends  toward  zero  (t*0)  all  the  £het9 
curves  coincide  with  the  uniform  local-oscillator  (ul  “  0)  case,  going  to  zero 
efficiency  in  the  zero-diameter  limit.  When  the  photodetector  diameter  becomes 
very  large  ( y  --  )  all  the  thet9  curves  approach  asymptotes  set  by  ul>  viz. 

( 18 )  reduces  to 

£het3  "  “l^1  "  exp(-4/uL2)j2/2  .  (23) 

From  (23),  the  large-7  Ehet9  asymptotes  can  be  shown  to  approach  zero  wtien 
uL-u  and  their  maximum  value  occurs  near  ul  “  1*2.  Finally,  because 

of  Bessel  function  oscillations,  the  7  ■»-  asymptotes  are  not  achieved 
monoton ica 1 1 y . 

To  those  accustomed  to  one-way  propagation  (i.e.,  laser  communication  rather 
than  laser  radar)  heterodyne  mixing  efficiency  calculations1-3.  Fig.  2  may  at 
first  present  some  puzzling  features.  It  is  well  known  that  unity  heterodyne 
mixing  efficiency  results  in  the  communication  problem  if  the  detector  diameter  d 
greatly  exceeds  the  focal-plane  width  of  the  signal  intensity  pattern  and  the 
local  oscillator  field  pattern  exactly  matches  the  signal  field  pattern  on  the 
detector.  It  might  seem,  therefore,  that  in  the  radar  calculation  thet9  will 
achieve  its  maximum  value  (if  not  unity)  when  the  detector  is  large  and  the  local 
oscillator  beam  (2)  is  matched  to  the  transmitter  beam  (1),  viz.  B/xult  <  1 
with  ul  *  ut  =  4aT/D-  Because  of  our  no-truncation  assumption  at  the 
transmitter  we  have  ut  £  1  »  and  thus  (22)  applies,  giving  thet9  £  for 

this  "matched"  case.  Even  the  maximum  thet9  value  over  7  and  ul  £  1  is 
still  less  than  the  m’xi-.uin  thet9  value  over  7  and  ul»  as  the  latter  occurs 
at  u l  >  1  • 

Physically  there  is  no  conflict  between  Fig.  2  and  the  laser  communication 
dictum  of  matching  the  local  oscillator  to  the  signal  over  the  detector. 

Equation  (1)  is  the  signal  field  leaving  the  transmitter  and  Eq.  (2)  is  the  local 
oscillator  field  impinging  on  the  detector.  For  an  unresolved  glint  target,  the 
target-return  field  impinging  on  the  detector  will  be  a  Bessel  function 
circular-pupil  diffraction  pattern.  Indeed  the  integral  011  the  right  in  (18)  is 
the  overlap  integral  of  this  Bessel  function  pattern  and  the  Gaussian  local 
oscillator  pattern.  That  thet9  has  a  maximum  value  at  fixed  3ql  as  d  is 


varied  is  therefore  exactly  like  the  laser  communication  mixing  efficiency 
calculation  for  a  plane  signal  wave  incident  on  an  objective  lens  being  mixed  on 
a  focal  plane  detector  with  a  Gaussian  local-oscillator  field3 ,  3  *  . 

4.2  SPECKLE  TARGET  CALCULATIONS 

For  far-field  propagation  with  the  radar  focused  at  infinity,  the  basic 
speckle-target  results  take  the  form  given  in  Section  1.2  with 

tT(p)  *  exp[-(*uTx)2/Bj  .  (24) 


T  2 

I  dy  y  exp[-(tuLy/4)  }  F(x,y) 

cR(p  »  -  [  -  1  1  '  <“> 

/  dy  exp[-(»u.y/4>2]  J,(»y) 

0 

and 

£het*  *  W5  ('V2)2  /’  d*  *  «r‘*>  V*>  '  (26) 

where 

x  -  D  |  p  j/AX.  (27) 

measures  the  transverse  distance  in  the  target  plane  relative  to  the  nominal 
far-field  diffraction  scale  of  the  objective  lens. 


uT  -  4aT/D  <  1 


(28) 


measures  the  transmitter  beam  relative  to  the  objective  diameter  lu^  <  1  is  our 
no  transmitter-beam  truncation  assumption). 


F(x,y) 


/  de  (2**)"1  J, (»z>, 
0 


(28) 


for 

z  -  (x2  +  y2  -  2xy  cos#)1/2  ,  (30) 

and 

g.(D|?|/XL)  -  for  i  -  T,  R,  (31) 

are  the  transmitter  and  receiver  angular  patterns  normalized  to  the  l/D 
diffraction  scale. 

As  was  the  cane  in  Section  4.1,  it  is  worth  considering  the  special  cases  of 
the  uniform  local  oscillator  (8/xult  >>  D,  and  the  untruncated  local 
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14.  Cohen,  S.C.,  "Heterodyne  Detection:  Phase  Front  Alignment,  Beam  Spot  Size 
and  Detector  Uniformity,"  Appl.  Opt.  1_4,  19S3-19S9  (1975). 


oscillator  t8/*uLY  <  1  and  u[,  <  D-  In  the  former  instance  we  find 

«R ( o  )  -  *2[l  -  )  j~2  G(y>,  (32) 

ChetS  ■  ( 2 w uT ) 2  /  dy  y  exp [ -( t p^y y ) 2/2 !  G(y),  (33) 

where 

y  -  |J|f/dL  (34) 

measures  transverse  distance  in  the  target  plane  relative  to  the  geometric 
projection  of  the  detector  diameter  onto  the  target  plane,  and 


Y  2 

G(y>  -  [/  dx  J^.x)  Ju (  2» xy )  j 


In  the  latter  case  we  obtain 


C R ( o  )  “  exp[-(*uLx)2/8]  , 


and 

‘».<Y  '  “T2“t.2/2I»T!  *  »t,2’  '  1371 

where  x  is  defined  in  (27). 

Let  us  first  use  the  foregoing  results  to  examine  the  radar's  spatial 
resolution.  The  composite  beam  pattern  ?t(p)Cr(p)  is  the  product  of  a  simple 
Gaussian  zt  with  e-2  radius  x  *  4/xuj  (equivalent  to  |pj  »  IL/xaj)  and  a 
complicated  cr  from  Eq.  (25).  In  Figs.  3  and  4  we  have  plotted  the  angular 
pattern  g^fx)  vs.  x  for  pi  *  0.5  and  1.7,  respectively,  and  a  variety  of  r 
values.  From  these  figures  the  following  characteristics  may  be  discerned.  For 
small  ul  and  large  y  the  beam  pattern  broadens,  approaching  the  Gaussian  (36) 
with  e-2  radius  x  »  4/*ul  (equivalent  to  |p|  •  XL/ra^  »  Lao[/f),  as 
expected.  For  ui  large  or  y  small,  the  beam  pattern  approaches  Airy  disk 
behavior,  with  a  deep  null  near  x’ »  1.22  (equivalent  to  jp|  ”  1.22  1L/D). 
Physically,  the  second  case  corresponds  to  a  uniform  local  oscillator  over  a 
detector  that  is  unresolved  by  the  objective  lens,  from  which  we  derive,  via  (32) 
with  y  <<  1, 

CR(p)  ♦  [2J1(»x)/xx] 2  .  (38) 

The  implications  of  the  above  results  are  really  not  surprising.  When  the 
local  oscillator  is  essentially  untruncated  Cr(p)  is  broader  than  the  iL/D 
objective-lens  diffraction  width.  When  the  uniform  local  oscillator  assumption 
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t  is  valid,  «r(p)  is  roughly  XL/D  in  radius.  For  tha  untruncated  case  we  obtain 

;  the  closed-form  expression 

« 

|  rres  -  XL/*(aT2  +  aL2)1/2  ,  (39) 

)  for  the  e"2  radius  of  cx<p>Cr<p)*  Generally  speaking,  increasing  aL  so 

i 

'  that  local  oscillator  truncation  is  incurred  improves  spatial  resolution 

somewhat,  but  the  resolution  is  never  better  than  about  0.7  XL/D  when  the  maximum 
ax  ■  D/4  is  employed.  (Note  that  the  reason  resolution  exceeds  the  1.22  XL/D 
•  Rayleigh  criterion  is  because  we  are  concerned  with  the  e“2  radius  rather  than 

J  the  radius  to  the  first  zero  of  the  composite  antenna  pattern). 

.  Let  us  turn  now  to  the  behavior  of  ehet3*  *hich  we  have  plotted  vs.  y  in 

Figs.  5  and  6  for  ux  *  and  !*°*  respectively,  and  various  ml  values. 

These  figures  share  many  of  the  features  found  in  the  ehet9  curves  from  Fig. 

2.  In  particular:  ehet3  has  a  maximum  value  for  each  ux*  ML  Paxr  as  T 
j  is  varied;  as  the  detector  diameter  tands  toward  zero  (y  *  0),  all  the  ehets 

j  curves  coincide  (corresponding  to  tne  uniform  local  oscillator  limit)  and 

» 

•  approach  zero  efficiency;  and  as  the  detector  diameter  becomes  very  large  (y  ♦  - ) 

|  all  the  Chet9  curves  approach  asymptotes  set  by  their  ml*  pt  values, 

|  which  for  fixed  ux  5°  to  zero  if  ml  *  11  or  PL  *  "•  The  speckle  target 

* 

'  mixing  efficiencies  are  appreciably  lower  than  the  glint-target  efficiencies, 

with  chets  “  0*28  being  the  maximum  attainable  value  within  the  parameter  sat 
we  have  plotted.  Part  of  this  reduction  is  forced  upon  us  by  the  no-truncation 
1  assumption  ux  £.  *3  a  simple  Fourier  optics  calculation  to  show  from  (8), 

(15),  (16)  that 

ehetS  <  4<XL)2  eaL  |eT'(0)|2/iD2  -  mt2/2  <  0.5  .  (40) 

Thus,  ehet3  *  0.28  actually  corresponds  to  a  56%  efficiency  vis  a  vis  the 
transmitter  limit  given  in  (40).  Note  that,  as  was  found  for  ehet9  in  Fi9* 

2,  this  maximum  is  not  reached  with  "matched"  operation  (i.e.,  8/xult  £  1  with 
pl  *  ux  £  D*  the  reason  being  that  the  preceding  conditions  do  not  best 
mat^n  the  local  oscillator  to  the  targec-return  signal  on  the  photodetector.  We 
ascribe  the  observed  ehet3  <  ehet9  behavior  to  the  partial  coherence  over 
the  objective  lens  of  the  target-return  field  from  the  resolved  speckle  target, 
as  opposed  to  the  complete  coherence  of  the  glint-target  return  over  this  lens. 
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5.0  NEAR-FIELD  CALCULATIONS 


It  was  pointed  out  quite  some  time  ago15  that  a  monostatic  coherent  laser 
radar  that  is  focused  at  infinity  will  suffer  a  very  significant  loss  in 
heterodyne  efficiency  when  viewing  targets  located  in  the  radar's  near  field. 
Gauss ian-beam  radar  calculations  show15  that  best  near-field  mixing  efficiency 
results  when  the  transmitter  beam  and  the  backpropagated  local-oscillator  beam 
have  their  geometric  optics  foci  at  the  target  plane.  Because  future  high- 
resolution  coherent  laser  radars  may  have  near  fields  extending  2  to  4  km,  the 
preceding  focus  effects  must  be  considered.  Owing  to  the  complexity  of  near¬ 
field  patterns  for  backpropagated  local  oscillator  beams  with  truncation,  we 
shall  limit  our  treatment  of  the  near-field  problem  to  the  untruncated  Maoi/d 
£  1,  4sl/D  <  1)  Gaussian  beam  regime. 

5.1  NEAR-FIELD  RESOLUTION 

The  near-field  propagation  untruncated  Gaussian  beam  patterns  are 

CT(p)  -  exp[-(*uTx)2/8(l  +2t2)]  ,  (41) 

and 

CR(p>  -  exp[-(«uLx)2/8(l  +  nR2) ]  ,  (42) 

where  ux»  ul»  ar.d  x  are  33  in  Section  4.0,  and 

[)T  -  kaT2  |R_1  -  L_1|/2  ,  (43) 


nR  »  kaL2  R-1  -  L_1  ;i 


(44) 


are  dimensionless  near-field  focusing  factors.  When  the  transmitter  and 
backpropagated  local  oscillator  have  their  geometric  foci  at  the  target  range 
(R  *  L),  Eqs.  (41)  and  (42)  reduce  to  the  far-field  formulas  (24)  and  (36), 
respectively,  as  expected.  The  near-field  spatial  resolution,  given  by  the  e"^ 
radius  of  Ct<p)Cr(p)  is  easily  found  to  be 


rres  “  ^A[*t2/(1  +  nT2>  +  aL2/(l  +  aR2)]  1/2 
which  simplifies  to  become  (39)  when  R  »  L. 


(45) 


15.  Sonnenschein,  C.M.,  and  Horrigan,  F.A.,  "S ignal-to-Noise  Relationships  for 
Coaxial  Systems  that  Heterodyne  Backscatter  from  the  Atmosphere,"  Appl. 
Opt.  10,  1600-1604  (1971). 


Because  2T  and  JJR  are  positive  for  R  *  L,  we  see  that  resolution 
degrades  in  the  near  field  when  the  geometric-optics  beam  foci  do  not  coincide 
with  the  target  range.  We  can  use  (45)  to  quantitatively  assess  whether  the 
foregoing  degradation  is  sufficient  that  high-resolution  near-field  radars  must 
use  dynamic  focusing  to  match  R  to  L  for  each  target  in  the  scanned  field  of 
view.  In  Fig.  7,  fpr  example,  we  plot  rres  vs.  range  L  for  X  *  10.6  um  with 
aT  m  aL  *  7.23  cm  and  various  choices  of  R.  These  *x,  a^  values  give  33 
urad  far-field  angular  resolution  when  collimated  beams  (R  ■  -)  are  employed,  in 
which  case  the  near  field  to  far  field  transition  begins  at  about  L  ■  1.5  km. 
Three  R  choices  are  shown  in  Fig.  7:  R  ■  2  km  (fixed  geometric-optics  focus  at  2 
km),  R  *  »  ( collimated-beam  operation),  and  R  *  L  (dynamic  focusing).  Over  a 
nominal  1  km  to  5  km  range  interval.  Fig.  7  indicates  that  fairly  comparable 
resolution  is  achie  'ed  by  the  three  R  choices  plotted,  implying  that  dynamic 
focusing  is  not  critical  for  good  resolution  in  this  case. 

5 .  ..  NEAR-FIELD  MIXING  EFFICIENCIES 

The  preceding  near-field  propagation  untruncated  Gaussian  beam  patterns  give 
rise  to  mixing  efficiencies 

ehet9  *  wL  /2(1  +  aR  )  ,  (46) 

for  a  glint  target,  and 

chet3  -  1/2[<1  +  nT2)/„T2  +  (l  ♦  aR2)/iiL2]  ,  (47) 

for  a  speckle  target.  Recall  that  we  are  referencing  our  heterodyne-mixing  areas 
to  the  objective  lens  area  *D2/4  to  obtain  these  efficiencies  (see  Eqs.  (9), 

(10),  (15),  (16,),  30  that  the  no-truncation  assumptions  uj  £  1  and  ul  $.  1 
imply  ehet^*  1  0.5  and  ehet3  i.  °*25.  These  maximum  Gaussian-beam 
efficiencies  are  achieved  when  ut  *  WL  *  and  R  *  L,  and  the  glint-target 
result  is  significantly  lower  than  the  best  far-field  truncated  local  oscillator 
ehet^  from  Fig.  2. 

As  was  found  for  near-field  resolution,  a?  and  aR  being  positive  for  R  * 

L  will  degrade  mixing  efficiency.  In  Fig.  8  we  plot  ef,et3  and  Ehets  vs* 

range  L  for  the  radar  considered  in  Section  5.1,  namely  X  ■  10.6  urn,  aT  *  aL 

*  7.23  cm,  with  R  *  2  km,  ■,  or  L,  and  ut  “  ^L  *  1*  (Because  a-j  *  a 

(46)  and  (47)  yield  e^et3  ”  ehet3/2*>  Whereas  the  spatial  resolution  was 

not  severely  affected  by  the  choice  of  geometric-optics  focal  distance  R,  Fig.  8 

shows  that  mixing  efficiency  is  a  strong  function  of  R.  Figure  8  shows  that 
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choosing  R  ■  2  km  achieves  better  than  80%  of  the  optimum  (R  *  L)  mixing 
efficiency  over  the  range  interval  from  1.2  km  to  5.5  km,  and  outperforms 
coll  imated-beam  (R  »  •• )  operation  over  the  range  interval  from  0.5  km  to  4  km. 
Thus,  dynamic  focusing  does  not  appear  to  be  necessary  here. 

6.0  CONCLUSIONS 

We  have  tried  to  provide  a  rigorous  construct  for  assessing  the  spatial 
resolution  and  heterodyne  mixing  efficiency  of  a  monostatic  coherent  laser 
radar.  The  basic  results  in  Section  3.0  have  rather  general  applicability,  e.g., 
they  could  be  used  to  study  performance  for  Cassegrain  telescopes  by  changing 
'  and  appropriately16.  The  results  of  Section  4.0  indicate  that 

“matched"  operation  of  an  untruncated  Gaussian  transmitter  beam  and  an 
untruncated  Gaussian  local  oscillator  beam  does  not  optimize  system 
performance17.  The  results  of  section  5.0  imply  that  fixed  focus  near-field 
operation  may  be  the  logical  approach  to  high  resolution  radars,  as  dynamic 
focusing  may  not  offer  performance  benefits  sufficient  to  warrant  the  concomitant 
system  complexity.  Finally,  we  should  point  out  that  atmospheric  turbulence 
effects  should  be  included  for  a  more  complete  treatment  of  spatial  resolution 
and  mixing  efficiency  of  large-aperture  radars. 


16.  Greene,  B.A. ,  Masters,  D.  ,  Rye,  B.J.,  and  Thomas,  E.L.,  "Antenna 

Considerations  in  Atmospheric  Backscatter  Coherent  Lidar,"  Technical  Digest 
of  Topical  Meeting  on  Coherent  Laser  Radar  for  Atmospheric  Sensing  (upt. 
Soc.  Am.,  Aspen,  1980)  paper  Tu03. 


17. 


Wang,  J.'i., 
( 1981  ). 


Heterodyne  Laser  Radar  Performance,"  Proc.  SPIE  300 ,  110-117 
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FIGURE  1.  MONOSTATIC  COHERENT  LASER  RADAR 
SYSTEM  STRUCTURE. 


FIGURE  2.  FAR-FIELD  GLINT-TARGET  MIXING  EFFICIENCY 
ehet9  vs-  NORMALIZED  DETECTOR  DIAMETER 
Y  FOR  VARIOUS  VALUES  OF  NORMALIZED  LOCAL- 
OSCILLATOR  INTENSITY  RADIUS  UL. 


0R  <K) 


FIGURE  3.  FAR-FIELD  RECEIVER  BEAM  PATTERN  gR(x)  vs 
NORMALIZED  TRANSVERSE  DISTANCE  x;  UL=*0.5 
AND  VARIOUS  Y  VALUES. 
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FIGURE  4.  FAR-FIELD  RECEIVER  BEAM  PATTERN  gR(x) 
vs.  NORMALIZED  TRANSVERSE  DISTANCE  x; 
Ut-1.7,  AND  VARIOUS  y  VALUES. 


FIGURE  5.  FAR-FIELD  SPECKLE -TARGET  MIXING  EFFICIENCY 
ehetS  vs>  NORMALIZED  DETECTOR  DIAMETER  y 
FOR  VARIOUS  VALUES  OF  NORMALIZED  LOCAL- 
OSCILLATOR  INTENSITY  RADIUS  PL;  uT=»0.7. 
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FIGURE  6.  FAR-FIELD  SPECKLE-TARGET  MIXING  EFFICIENCY 
ehetS  vs.  NORMALIZED  DETECTOR  DIAMETER  y 
FOR  VARIOUS  VALUES  OF  NORMALIZED  LOCAL- 
OSCILLATOR  INTENSITY  RADIUS  UL;  u^l .  0  . 
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Abstract 


Coherent  laser  radars  used  for  3-D  Imaging  and  aerosol/pollutant  sensing  collect  target 
returns  which  comprise  a  range-spread  speckle  process.  Thts  paper  addresses  a  class  of 
waveform  deslgn/waveforra  evaluation  problems  for  such  radars.  The  maximum-likelihood  (ML; 
processors  and  Cramer-Rao  (CR)  performance  bounds  are  developed  for  short  pulse,  chirped 
pulse,  sinusoidal  amplitude-modulated  pulse,  and  sinusoidal  frequency-modulated  pulse  wave¬ 
forms.  All  four  of  these  waveforms  are  shown  to  have  nearly  identical  ultimate  range 
accuracies.  However,  range  accuracy  is  not  the  only  performance  criterion  of  Interest; 
range  resolution,  ambiguity,  and  anomaly  must  also  be  considered.  Analysis  and  discussion 
of  some  of  these  aspects  of  waveform  selection  are  Included  In  the  paper. 

Introduction 


Coherent  laser  radars  are  now  being  developed  for  a  variety  of  Imaging,  detection,  and 
remote  sensing  applications.  In  two  Important  application  areas,  3-D  imaging  and 
aerosol/pollutant  sensing,  t.ie  target  returns  comprise  a  range-spread  speckle  process. 

This  paper  addresses  fundamental  performance  limitations,  and  specific  performance  results 
for  targets  of  this  class.  We  begin  by  examining  the  range  accuracy  for  isolated  speckle 
targets  achievable  with  four  prototypical  pulse  waveforms  and  maximum-likelihood  (ML)  pro¬ 
cessors.  Cramer-Rao  (CR)  performance  bounds  are  developed  which  show  that  all  four  of 
these  waveforms  have  nearly  Identical  ultimate  range  accuracies.  However,  In  range-spread 
target  environments  range  accuracy  Is  not  the  only  performance  criterion  of  interest. 

Range  resolution,  which  relates  to  the  Instantaneous  range  swath  sampled  by  the  pulse  at  a 
given  time  Instant,  Is  also  important.  In  addition,  range  ambiguity,  arising  from  the 
structure  of  the  waveform,  and  range  anomaly,  created  by  noise,  must  also  be  considered. 
Some  general  features  of  these  performance  measures,  as  well  as  some  specific  system  calcu¬ 
lations  will  be  presented. 

Radar  signal  model 

The  model  we  employ  In  our  analysis  Is  the  range-spread  generalization  of  the  model 
used  In  our  previous  work.1*2  We  assume  a  compact  coherent  laser  radar  whose  operation  Is 
unaffected  by  atmospheric  turbulence  and  whose  Intermediate-frequency  (IP)  received  signal 
has  complex  envelope 


_r(  t )  ■  £(  t )  +  n  ( t ) , 


(1) 


where  Is  $he  target  return  and  ji(t)  Is  local-oscillator  shot  noise.  The  target 

return  for  a  single  pulse  transmission  obeys 

m 

X(t)  •  E1/2  /  dz  s.(t  -  2z/c)  /dp"'  Tip1  ,z,t-z/c)  £2(p~,z),  (2) 

0  ~ 

where:  E  Is  the  transmitted  pulse  energy;  s_(t)  Is  the  normalized  complex  envelope  of  the 
transmitted  pulse;  T(p,'z,t)  Is  the  random  speckle-target  field  reflection  coefficient  at 
transverse  coordinate  JT*  range  z  and  time  t;  and  ^(p’,z)  Is  the  common  normalized  spatial 
beam  pattern  of  transmitter  and  back-propagated  local  oscillator  at  range  z.  The  shot 
noise  jn(t)  Is  a  zero-mean  clrculo-complex  white  Gaussian  noise  with  spectral  density  hvo/n 
In  terms  of  the  photon  energy  and  the  detector  quantum  efficiency. 

We  shall  assume  that  TlT'.z.t)  Is  a  c 
function 

<T(7i'  ,zi,  ti)  T(p2'  , 

.  S(z\-Z2) 


Ircuio-complex  Gaussian  process  with  covariance 

z2>  z2'*  >  *  T(2l)  Svpi'-pj' ) 

exp[-( ti-to) 2  /  tc2]. 


(3) 
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Physically,  this  model  corresponds  to  a  stationary  spatially-resolved  speckle  target  with 
range-dependent  average  reflection  strength  T(z)  and  coherence  time  tc.  For  the  par¬ 
ticular  case  of  an  aerosol  target  we  Identify 

3a(t,z)  ■  *  T(z)  (4) 

as  the  atmospheric  backscatter  coefficient  (m-1  sr-1  units)  at  range  z,  and  tQ  as  the 
target-return  correlation  time  arising  from  the  aerosol's  velocity  fluctuations. 

Por  a  stationary  range-spread  hard  target  we  can  usually  assume  tc  *  -  and  treat  T(z) 
as  the  target  characteristic  of  Interest.  Note  that  the  isolated  range-L  stationary 
hard-target  model  used  in  Refs.  1  and  2  corresponds  to  Eq.(3)  with  tc  *  *  and 

T ( z )  -  (e/»)  d(z-L)  (5) 

for  p  the  target's  diffuse  reflectivity. 

Range  accuracy 

Suppose  the  radar  Is  used  against  the  Isolated  range  speckle  target  Just  described,  and 
the  IF  signal  for  a  single  pulse  is  to  be  processed  to  yield  an  estimate  of  the  range  L. 

The  maximum-likelihood  range  estimate,  l,  Is  the  z  value  that  maximizes^ 

£(  z)  ■  |  /  dt  jr(t)  s_(  t-2z/c)*  j  2.  (6) 

Physically,  this  processor  can  be  realized  with  a  matched-filter/envelope-detector/peak 
detector  cascade. 

A  range  estimate  L  Is  said  to  be  unbiased  if  Its  average  value  (over  the  shot  noise  and 
target  speckle  ensembles)  <L>  equals  the  true  range  L,  no  matter  what  value  L  was.  All 
unbiased  estimates  have  mean-squared  estimation  errors  that  obey  the  Cramer-Rao  Inequality3 

d£  =  <( L  -  L)2  >  • 

5  2 
c2  (CNR  +  1)  at 


8  CNR2  [4  4- 

In  Eq.  (7):  c  Is  the  speed  of  light; 

CNR  ■  nEo  Ar  e-2^  /  hv0  *L2 
Is  the  carrier-to-nol3e  ratio,  where  a  Is  the  atmospheric  extinction  coefficient  and 

Ar  -  UL)2  e2S“  /dp~  |  _£(  ~,  L)  |  **  (9) 

is  the  effective  heteodyne-receptlon  area; 
and 


2 

°t  * 

/dt 

t2  |  s_(  t) 

|  2  _  [/dt  t  |  s_(t)  )  2]2> 

(10a) 

2 

°u  * 

/dt 

|  s.(t)  |  2 

+  [/dt  _s(t)  ^(t)*]2, 

(10b) 

*tw  ■ 

Im 

[/dt  _a(t) 

• 

s(t)*]  +  J/dt  t  |  s.(t)  (  2  /dt  s(t)3.(t)*, 

(10c) 

are  parameters  determined  by  the  waveform  s_(t)  that  the  radar  employs.  The  parameters  <jt 
and  ou  are  the  effective  root-mean-square  time  duration  and  radian-f requency  bandwidth  of 
_s(t),  respectively.  The  parameter  measures  the  degree  of  range-Doppler  coupling  in  the 
waveform  (see  below)  with  ■  0  corresponding  to  zero  coupling.  Note  that  when  *  0 
and  CNR  >>  1  Eq.  (7)  reduces  to 

■5L  >  c/23/2  Clsal/2  Xii,  (ii) 


(7) 


(8) 
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l.e.,  the  ultimate  range  accuracy  goes  Inversely  with  signal  bandwidth  and  the  square  root 
of  the  carrier-to-nolae  ratio. 

Let  us  examine  the  waveform  dependent  factor  In  the  CR  bound,  viz. 

S  ■  /  C  at  ®w  "  xtm  31,/2  (12) 

for  four  prototypical  s(t):  1)  a  short  Gaussian  pulse,  2)  an  upchlrped/downchlrped 
QauJ'lan  pulse,  3)  a  sinusoidal  amplitude-modulated  (am)  Gaussian  pulse,  and  4)  a  sinu¬ 
soidal  frequency-modulated  (fm)  Gaussian  pulse.  The  specific  waveforms  are  as  follows. 

Short  Gaussian  pulse 

Here  we  assume 


s_(t)  *  K-i  exp  C-(fWt)2/4]  (13) 

where  Kj  is  a  normalization  constant,  W  is  the  Fourier  transform  bandwidth  (In  Hz)  and 
4/»W  is  the  time  duration  of  the  pulse*. 


(Jpchlrped/downchlrped  Gaussian  pulse 


In  this  case  we  assume 

aU)  -  K2  exp  [-4t2  /  T2  -  J(2x)l/2  Wt] 
•  exp  [J2*Wt  j  t  |  /T] , 


(14) 


where  T  Is  the  time  duration,  W  is  the  chirp  bandwidth,  and  WT  >>  1.  Mote  that  we  have 
chosen  the  upchlrped/downchlrped  waveform  to  avoid  the  very  strong  range-Doppler  coupling  In 
the  ordinary  chirped  pulse.  The  latter  causes  significant  problems  for  non-statlonary 
targets  whose  range  and  Doppler  shift  must  be  estimated  from  a  single  pulse. 


Sinusoidal  am  pulse 

In  this  case  we  assume 

s.(t)  -  K3  exp  C-4t2  /  T2  ]  [1  +  mcos(*Wt)],  (15) 

where  T  Is  the  pulse  duration,  m  Is  the  modulation  strength,  W  Is  the  modulation  bandwidth, 

and  WT  >>  1. 

Sinusoidal  fm  pulse 

Here  we  assume 

S,(t)  •  Kn  exp  [-4t2  /  T2  +  J»(W/j3)cos(wst)  ],  (16) 

where  T  Is  the  pulse  duration,  W  is  the  full  width  of  instantaneous-frequency  swing,  a>3  Is 

the  sinusoidal  modulation  radian  frequency,  and  WT  >>  1,  «,T  »  1. 

S-factor  value 


All  four  of  the  preceding  waveforms  have  the  same  nominal  bandwidth  and  the  last 
three  have  the  same  nominal  time  duration  T.  Thus,  by  comparing  the  values  these  waveforms 
yield  for  S  In  Eq.  (12)  we  can  see  how  their  ultimate  range  accuracies,  l.e.,  their 
respective  CR  bounds,  compare.  We  find  that 


2/w 

-  0.64, 

short  pulse 

(j2/2  -  I)-!/2 

■  0.50, 

chirped  pulse 

(2  +  m2) l/2/mw 

*  0.55, 

for  m"l,  am  pulse 

2l'/2/» 

-  C.45, 

fm  pulse. 

(17) 


For  the  pure  Gaussian  cases,  bandwidth  and  time  duration  are  measured  full  width  to  e"2 

intensity  points. 
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(Our  choice  of  m*l  for  the  am  pulse  Is  dictated  by  the  desire  for  envelope  detector 
compatibility).  It  Is  apparent  from  these  numbers  that  all  four  waveforms  have  nearly 
Identical  range  accuracies.  Accordingly,  the  accuracy  achievable  with  ML  processing  cannot 
serve  as  a  reason  for  choosing  one  waveform  over  another. 

Range  resolution 


A  second  waveform  performance  criterion  of  Interest  la  Its  range  resolution  zres* 
Roughly  speaking,  zres  Is  the  fanS*  swath  which  can  give  rise  to  a  target  return  at  some 
particular  time,  e.g. ,  from  (2)  It  follows  that  zres  »  cT/2  for  a  rectangular  pulse  of 
duration  T.  Range  resolution  plays  a  role  In  distinguishing  discrete  targets  located  at 
various  ranges  from  the  radar.  For  range-spread  targets  such  as  aersols,  zres  I®  6^e 
awath  over  which  the  backscatter  profile  is  smeared.  Finally,  as  will  be  discussed  later, 
range  resolution  enters  Into  the  anomaly  calculation. 

To  Illustrate  the  behavior  of  zres  mentioned  above,  consider  the  case  of  aerosol 
sensing  using  the  chirped  Gaussian  pulse 

js ( t )  -  K  exp  [-4t2  /  T2  -  J*Wt2  /  T],  (IS) 

which  has  duration  T  and  chirp  bandwidth  W,  with  WT  >>  1  assumed.  We  shall  assume  far- 
fleld  propagation  using  a  collimated  Gausslan-beam  transmitter/local  oscillator  beam  of 
Initial  a-2  Intensity  radius  a.  We  shall  also  assume  matched-filter/envelope-detection  Is 
performed  for  range  z,  l.e.,  the  receiver  generates  f(z)  given  by  (6).  It  follows  that  the 
carrler-to-nolse  ratio  for  range  z  can  be  written  In  the  form 


where 


CNR(z)  •  /  dz'  W(z')  K(z-z'), 

0 

Q( z)  ■  [n  E  c  T  .3/2  a2  as(i,Z)/hv0  4z2]  exp[-2  /  %(z’)  dz'] 


(19) 

(20) 


Is  the  Infinite  resolution  aerosol-return  CNR  profile,  and 


2  2 
p(z)  *  g(*  Zres  /  8)“1/2  eXp  (_8  z2/ Zres) 


(21) 


is  the  effective  probe  function  of  the  chlrped-transmltter/matched-f liter  receiver. 


In  Eq.  (21)  the  range  resolution  Is  given  by 

cT/21/2 


ZRES  ‘ 


[1  +  (iWT/4)2  /  (1  +  ^/atc2)]1/2’ 
g  »  /dz  F(z)  *  [1  +  T2/4tc2  +  (tWT/4)2]-!/2. 


(22) 

(23) 


For  example,  suppose  T  ■  10us,  W  »  20MHz  and  t0  »  lus;  we  find  that  zRES  »  68.8ra.  Were 
there  no  aerosol  decorrelation,  l.e.,  tc  »  -,  these  same  T  and  W  values  yield  much  better 
resolution,  zres  “  13.5m.  On  the  other  hand,  for  a  smoothly  varying  aerosol  distribution. 


CNR(z)  -  Q(  z )  g 


(24) 


and  g  varies  by  less  than  0.1*  for  T  and  W  as  above  when  tc  changes  from  infinity  down  to 
lus.  Physically,  when  t0  <  -  the  loss  of  correlation  over  the  pulse  duration  T  spoils  the 
resolution  by  spoiling  the  pulse-compression  In  the  matched  filter.  Yet,  because  the 
broadened  probe  function  now  collects  aerosol  returns  from  a  larger  range  3wath,  there  Is 
almost  no  Immediate  I033  In  CNR. 


Ranp;e  ambiguity 

The  discussion  thus  far  of  range  accuracy  and  range  resolution  presumes  that  in  any 
problem  there  Is  a  single  range  Interval  of  Interest.  In  particular,  in  the  accuracy  work 
we  were  Interested  In  an  Interval  of  width  25t  about  L,  and  In  the  resolution  work  we  were 
Interested  In  an  Interval  of  width  zres  about  z.  Many  radar  waveforms  have  range  ambi¬ 
guities,  l.e.,  for  a  target  at  range  1  there  are  a  variety  of  z  values  which  yield  nearly 
Identical  noise-free  matched-filter  aspects  outputs. 

Iq(z)  «  |  /  dt  _s(  t-2L/c )  s_  (t-2z/o)  *  j  2. 
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(25) 


Because  of  noise  effects,  the  radar  trill  be  hard  pressed  to  correctly  sort  out  which  of 
these  ambiguous  range  values  to  associate  with  the  target.  In  this  regard  the  sinusoidal 
am  and  fm  pulses  we  considered  under  range  accuracy  are  at  a  disadvantage  trlth  respect  to 
the  short  pulse  and  chirped  pulse  waveforms  considered  therein.  For  the  sinusoidal  am 
pulse  we  have 

fo(L-ck/W)  -  exp(-l6  k2/W2  T2)  (26) 

for  k  an  Integer,  which  Implies,  because  WT  >>  1,  range  ambiguities  exist  every  c/W. 
Likewise,  for  the  sinusoidal  fa  pulse  and  k  an  Integer  we  obtain 

*  foC1*  -  *kc/wg)  *  exp  (-16*2  k2/*^2  T2) ,  (27) 

Implying  range  ambiguities  every  nc/wj,  as  wjT  >>  1. 

Range  anomaly 

Even  when  the  radar's  waveform  s_(t)  Is  not  prone  to  range  ambiguity,  a  single-pulse 
range  estimate  from  a  speckle  targe?  will  not  have  an  accuracy  well  approximated  by  the  CR 
bound  (7).  Indeed,  the  range  estimate  may  be  anomalous,  l.e..  It  may  lie  outside  of  the 
2 RES  Interval  that  contains  the  target.  This  Is  the  well  known  threshold  problem  In 
nonlinear  estimation  theory14.  Physically,  a  single-pulse  speckle  target  return  may,  owing 
to  destructive  Interference,  be  weak  enough  that  the  ML  processor  chooses  a  range  estimate 
corresponding  to  a  noise  peak  that  may  be  arbitrarily  far  removed  from  the  correct  range- 
resolution  Interval.  Por  an  a  priori  target  range  uncertainty  region  of  width 

42  ■  M  2RE5  (28) 

no  Doppler  uncertainty,  and  a  radar  waveform  without  ambiguity  problems  the  probability  of 
range  anomaly  on  a  single  pulse  is  approximately^ 

Pr( A)  -  (In  M  -  (2M)“1  ♦  0.577]  /CNR.  (29) 

Because  the  root-mean-square  range  estimation  error  given  a  range  anomaly  has  occurred  Is 
roughly  A2/61/2,  the  single-pulse  range  accuracy  Is  actually  dominated  by  these  events, 
l.e.,  Pr(A)  az/61'2  typically  exceeds  the  CR  bound  (7). 

In  order  to  circumvent  the  anomaly  problem,  a  sequence  of  N  Independent  pulses  can  be 
employed.  The  anomaly  probability  now  becomes 

Pr(A)  •  M  (l+CNR)N/(»N)l/2  CNR  ( l+CNR/2) 2N-1 ;  (30) 

It  follows  for  reasonable  N  and  CNR  values  that  N-pulse  range  accuracy  can  approach  the  N- 
pulse  CR  bound_ 

iL  •  c  3  (CNR  +  1)1/2  /  83/2  cNh  N 1/2 .  (31) 
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ABSTRACT 

Two  laser  radar  signal  models  are  developed,  one  for  a 
Doppler-shifted  target  return  with  randomly  time-dependent 
speckle  and  short-term  laser  frequency  Instability,  and  one 
for  a  clutter  return  with  the  above  characteristics  augmented 
by  a  random  macro-motion  Doppler  shift.  Target-return 
velocity  resolution  for  a  Doppler  radar,  and  target-return 
range  resolution  for  a  chirped  radar  are  derived  from  the 
former  model.  C 1  utter- 1 i mi  ted  false  alarm  probabilities  are 
derived  from  th*  latcer.  Experimental  Doppl er-radar  returns 
from  a  moving  flame-sprayed  aluminum  calibration  plate  are 
analyzed  for  comparison  with  the  target  speckle  model. 

Experimental  Oopp 1 er-radar  tree  clutter  returns  are  analyzed 
for  comparison  with  the  clutter  model. 

1.0  INTRODUCTION 

Velocity-imaging  Doppler  laser  radars,  and  range-imaging  laser  radars  which 
use  high  time-bandwidth  (TW)  product  waveforms,  both  employ  pixel  dwell  times  far 
in  excess  of  the  pulse  durations  of  pulsed-lmager  laser  radars.  The  long  dwell 
times  make  the  performance  of  the  former  systems  subject  to  degradations  arising 
from  laser  frequency  instability  and  target  decorrelation  to  which  the  latter 
systems  are  largely  Immune.  These  degradations  must  be  Included,  along  with 
speckle  effects  and,  for  Doppler  radars,  clutter,  in  analyses  of  resolution  and 
signal -to'-noise  ratio  for  such  systems.  In  this  paper,  preliminary  models  for 
the  preceding  effects  and  resulting  performance  analyses  are  presented.  Data 
'.ollected  with  the  MIT  Lincoln  Laboratory  Doppler  radar  are  used  to  test  some  of 
the  speckle  and  clutter  characteristics  assumed  In  the  models. 

Section  2.0  presents  the  radar  system  structure  to  be  analyzed  and  develops 
two  radar  signal  models.  The  first  model  describes  a  Doppler-shifted  target  re¬ 
turn  with  randomly  time-dependent  speckle  and  short-term  laser  frequency 


This  work  was  sunoorted  by  the  Department  of  the  Air  Force  Contract  F1D628-30-C-0002 
and  the  II. 5.  Amy  Research  Office  Contract  "AAORp-.r i-<-no?2.  "The  U.F,  C-overnment  assumes 
no  responsibility  for  the  information  Presented. ’’ 
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instaOT  1  it>.  The  second  model  characterizes  a  clutter  return  witn  tne  above 
features  augmented  by  a  random  macro-motion  Doppler  shift.  Performance  analyses 
using  these  models  are  given  in  Section  3.0.  Velocity  resolution  for  a  Doppler 
radar,  and  range  resolution  for  a  chirped  radar  are  derived  from  the  target 
return  model.  Cl  utter- 1 1  mi  ted  false  alarm  probabilities  are  derived  from  the 
clutter  model.  Section  4.0  uses  experimental  results  to  examine  the  validity  of 
the  models.  Spatially  averaged  intensity  histograms  of  experimental 
Doppler-radar  returns  from  a  moving  f 1 ame-sprayed  aluminum  calibration  plate  are 
analyzed  for  comparison  with  accepted  speckle  behavior.  Velocity  histograms  of 
experimental  tree  clutter  are  analyzed  for  comparison  with  the  assumed  clutter 
model  . 

2.0  SYSTEM  MODELS 

The  basic  system  structure  we  shall  assume  Is  shown  schematically  In 
Fig.  1.  It  Is  a  monostatic  shared-optics  coherent  laser  radar  using  a  single 
laser  for  both  tne  transmitter  and  local  oscillator  functions.  Long  dwell  time 
operation  without  transmitter  modulation  is  used  for  vel ocl ty-1 magi ng  operation. 
Long  dwell  time  operation  with  a  chirped  (Internally  or  externally  modulated) 
transmitter  Is  used  for  range-imaging  operation.  In  either  case,  an  externally 
generated  frequency  offset  Is  used  to  obtain  the  local  oscillator  beam  for 
unaliased  heterodyne  detection  of  the  radar  return. 

As  in  our  earlier  studl  es 1 ■ 2> 3» “  the  fundamental  quantity  of 
interest  Is  the  complex  envelope  of  the  Intermediate-frequency  (IF)  received 
signal  for  a  single  pixel  dwell  time,  given  by 

£(t)  *  jv(t)  +  njt)  (1) 
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2.  Papurt,  D.  M, ,  Shapiro,  J.  H.,  and  Harney,  R.  C. ,  "Atmospheric  Propagation 

Effects  on  Coherent  Laser  Radars,"  Proc.  SPIE  300.  86-99  (1981). 

3.  Papurt,  0.  M. ,  Shapiro,  J.  H. ,  and  Lau,  S.  T. ,  "Measured  Turbulence  and 

Speckle  Effects  In  Laser  Radar  Target  Returns,"  Proc.  SPIE  415,  In  press 
(1983).  - 

4.  Shapiro,  J.  H.,  "Resolution  and  Accuracy,  Ambiguity  and  Anomaly  for 
Range-Spread  Speckle  Targets,"  Proc.  SPIE  415.  In  press  (1983). 


in  terms  of  tne  radar  return  t )  ana  the  local  oscillator  snut  noise  Mt).  The 
excellent  spatial  resolution  of  laser  radar  equipment  makes  it  reasonable  to 
assume  tnat  /(t)  is  either  a  target  return  or  a  clutter  return,  but  not  a  combin¬ 
ation  of  tne  two.5  The  1 ocal -osci 1 1 ator  shot  noise  M t )  is  a  zero-mean 
ci rculo-compl ex  white  Gaussian  noise  process  of  spectral  height  hv0/n  for  hv0  the 
photon  energy  at  the  laser  wavelength  and  n  the  quantum  efficiency  of  the 
photodetector . 

2.1  TARGET-RETURN  MODEL 

We  assume  a  compact  conerent  laser  radar  whose  operation  is  unaffected  by 
atmospneric  turbulence,  but  is  subject  to  atmospheric  extinction  with  coefficient 
a  m'1.  A  range-spread.  Doppler-shifted,  time-dependent  speckle  target  will  then 
give  rise  to  a  single  pixel  target  return  complex  envelope1* 

l( t)  *  E1/2  exp { j2 if Qt )  /  dz  s(t-2z/c)  exp[-j( ♦( t-Zz/c )- ♦( t ) ) ] 

0 

./  d7  l(7,z,t-z/c)  e'aZ  ^2(7,Z).  (2) 

In  Eq.  (2):  E  is  the  energy  transmitted  during  the  dwell  time;  fg  *  2vz/X  is 
the  Doppler  frequency  shift  associated  with  the  nominal  longitudinal  target 
velocity  vz  at  the  radar  wavelength  A;  £(t)  is  the  normalized  complex  envelope 
of  the  transmitted  field  neglecting  frequency  instability  effects;  ^ ( t )  is  the 
random  phase  shift  of  the  transmitted  field  that  is  due  to  laser  frequency 
instability;  _T( o ,  z,  t)  is  the  random  speckle-target  field  reflection  coefficient 
(neglecting  the  net  Doppler  shift  fg)  at  transverse  coordinate  o  *  (x,y)  range 
z  and  time  t;  and  _C(C,z)  is  the  common  normalized  free-space  spatial  beam  of  the 
transmitter  and  local  oscillator  at  range  z.6 

For  velocity  imaging  we  shall  take  £(t)  to  be  a  transform  limited  Gaussian 

pulse 

s(t)  •  (8/iT2)1M  exp  (  -4t  2/T  2) ,  (3) 


5.  Harney,  R.  C.,  “Design  Considerations  for  the  Infrared  Airborne  Radar  (IRAR) 
MTI  Subsystem,"  Project  Report  TST-26,  Lincoln  Laboratory,  M.I.T,  (July 
198U  ) . 

6.  The  assumption  of  common  beam  patterns  has  been  made  for  simplicity.  An 
examination  of  transmitter  and  1 ocal -osc i 1 1 ator  beam  pattern  effects  on  the 
spatial  resolution  and  heterodyne  mixing  efficiency  of  a  monostatic  radar  is 
given  in  Shapiro,  J.  H.,  Dardzinski,  V.  E.,  and  Tung,  E.  W.,  "Coherent  Laser 
Radar  Antenna  Patterns  and  Mixing  Efficiencies,"  Active  Systems  Specialty 
Group,  IRIS  (November  1983). 
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wnose  dwe 1 1  time  (full  widen  to  e"*  points  of  ( t  >  J  ^ )  T  Is  tans  of  usee  to’ 
provide  satisfactory  Doppler  resolution.  In  this  case  it  is  reasonable  to  assume 
that  tne  target  is  unresolved  In  range,  viz.  T(J.z-t)  in  Eq.  (2)  is  a  zero-mean 
c i rcul o-compl ex  Gaussian  process  witn  covariance  function 

<X(  P  1  «2  1  •  £  1  )i  (<52>22>^2)>  *  (*2  P/») 

•S(pi-P2)  4(Zl*t2)  4(*1*L) 

•exp[-(t  1-t 2) 2/tc2].  (4) 

Here  0  is  the  (assumed  uniform)  diffuse  ref lectivit”  of  the  target,  l  is  the 
nominal  target  range,  and  tc  Is  the  coherence  time  of  the  target’s  speckle 
process . 

For  range  imaging  we  shall  assume  fQ  •  0  to  avoid  range/Doppler 
ambiguity7,  and  take  $.(t)  to  be  the  chirped  Gaussian  pulse 

s(t)  -  (8/ »T  z) 1/4  exp[-4t  Z/T  2  -  j*Wtz/T]  (5) 

of  duration  T  and  chirp  bandwidth  W,  with  TW>>1.  In  this  case  the  target  may  be 
resolved  in  range,  so  we  modify  (4)  to  become 


<K0l.2l.ti)  ]_  (  P2»2  2,t  2)^  •  AZT(Zi) 

•  4  ("0 1  -"0  2 )  4(  z  i*z  2)  exp[-  (t  i*t  2)  2/tc  2],  (6) 

in  terms  of  the  range-dependent  average  reflection  strength  T(z).  For  the 
Important  example  of  an  aerosol  target  we  identify  . 

8S( *,z)  •  *  T  (z)  ( 7 j 

as  the  atmospheric  backscatter  coefficient  (m'l  sr*1  units)  at  range  z. 


In  real  chirped  systems  the  ranye/Doppl er  ambiguity  can  be  resolved  by  using 
a  composite  waveform  consisting  of  an  up-chirp  segment  and  a  down-chirp 
segment,  see,  e.g.,  Skolnik,  M.  I.,  Introduction  to  Radar  Systei-.s 
(McGraw-Hill,  New  York,  1980),  pp.  81-92. 
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It  still  remains  for  us  to  characterize  the  frequency  ins  caD  i  1 1  ty  term 
exp  C- j ( ♦( t-2z/c )- t( t ) )  ]  appearing  In  Eq.  (2),  which  Is  an  effect  not  previously 
treated  in  our  pulsed  imager  work.1.2*3  Physically,  this  term  arises  from 
the  beat  between  the  range  delayed  random  phase  shift  of  the  transmitter  beam  and 
the  undelayed  random  phase  shift  from  the  same  laser  used  as  the  local 
oscillator.  Typical  CO 2  lasers  used  in  coherent  laser  radars  have  very  narrow 
short-term  llnewldths  with  slowly  (and  randomly)  varying  center  frequencies  that 
lead  to  appreciably  larger  long-term  llnewldths.  fo  model  this  behavior  we  shall 
write  the  random  phase  process  as 

t 

♦(t)  *  2*  /  f j( t)dt  (8) 


In  terms  of  an  Instantaneous  frequency  fj(t)  that  is  a  zero-mean  stationary 
Gaussian  random  process  with  covariance  function 

<fj (t+t)f j (t)>  *  of2  exp ( - |  r|/tf ).  (9) 

In  this  formulation,  of  Is  the  root-mean-square  (rms)  long-term  laser 
bandwidth,  and  l/*tf  gives  the  3d8  short-term  laser  bandwidth,  where 
oftf>>l  is  assumed. 

From  Eqs.  (8)  and  (9)  It  Is  a  straightforward  but  tedious  calculation  to 
compute  the  mean  function  and  correlation  function  of  the  frequency  Instability 
term  for  range  z 


w,(t)  s  exp(-j( ♦(t-2z/c)  -♦(t))].  (10) 

Of  particular  Interest  are  the  following  approximate  results  for  the  spectral 
dens  1 ty  of  w 2 ( t ) 
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Sww<n  ;  dr  <*2(fr)w*(t)>  e*j2,fT 

m  m 


2tf (4»ofz/c) 2 

(4„o  z/c)2]  5(f)  ♦  — - - T- -  , 

(2«ftf )2  ♦  1 


for  2afz/c<<l,  (11a) 


for  2z/ctf<<l,  2ofz/c>>l,.  (lib) 
(4.0f2)*1/2  exp  ( -  f  2/4  <jf  2 )  , 

for  2z/ctf>>l,  (11c) 


Mhlcn  have  been  plotted  In  Fig.  2.  Equation  (11a)  applies  wnen  the  2z/c  sec 
propagation  delay  Is  Insufficient  to  resolve  the  long-term  laser  llnewldth;  the 
resulting  spectrum  comprises  a  spike  plus  a  Lorentzlan  of  3dB  bandwidth  l/»tf. 
Equation  (lib)  corresponds  to  the  normal  operating  regime  for  laser  radars  In 
which  the  propagation  delay  resolves  the  long-term  llnewldth  but  not  the 
short-term  llnewldth;  here  the  spectrum  Is  Gaussian  with  rms  bandwidth 
2of ( z/ct f ) l/2 ,  which  Is  much  broader  than  l/»tf  but  much  narrower  than 
of.  Finally,  Eq.  (11c)  describes  very  long  path  operation  for  which  the 
propagation  delay  resolves  the  short-term  llnewldth;  here  the  spectrum  Is 
Gaussian  with  rms  bandwidth  because  the  delayed  transmitter  beam  has  a 

frequency  Instability  that  Is  statistically  Independent  of  that  of  the  local 
oscl 1 1 ator . 

2.2  CUJTTER-RETURN  MODEL 

An  Important  consideration  In  assessing  the  utility  of  a  Doppler  radar  as  a 
moving  target  Indication  (MTI)  system  Is  Its  ability  to  distinguish  a  slowly 
moving  hard  target  from  background  clutter  returns  caused  by  windblown  trees, 
etc.5  Equations  (1)  and  (2)  are  valid  for  background  clutter,  with  T(7,  z,  t) 
now  being  the  complex  field  reflection  coefficient  of  the  clutter  at  transverse 
coordinate  7,  range  z,  and  time  t.  For  velocity  Imaging  with  the  long  dwell  time 
transmitted  pulse  (3),  we  hypothesize  that  X(7,z,t)  will  be  such  that  a  clutter 
return  jr  ( t )  will  take  the  form 


1/2  ^ 

£(t)  «  £  c  ‘  s_(  t -2  L/C  )  ;  t  )£c  ( t )  exp  [  j  2  «  J  fc(t)dt],  (12) 

where  £c  is  the  average  received  clutter  energy,  s_(t)  is  the  normalised 
transmitted  field  complex  envelope  neglecting  frequency  instability  (given  by 
Eq.  (3)),  L  is  the  nominal  range,  and  W[_ ( t )  1s  the  frequency  instability  term 
(given  by  Eq.  (10)).  The  new  features  in  Eq.  (12)  are  a  clutter-speckle  process 
£c(t),  taken  to  be  a  zero-mean  uni ty-vari ance  stationary  cl rculo-compl ex 
Gaussian  random  process  with  spectral  density 

S  (f)  »  (2»o  2)*1/2  exp  ( -f  2/2  o  2).  (13) 

ycyc  c  c 

and  a  clutter  macro-motion  Instantaneous  Doppler  shift  fc(t),  taken  to  be  a 
Gaussian  random  process  that  is  statistically  Independent  of  £c(t)  with  mean 
mf  c ,  variance  <jfC2,  and  correlation  time  far  in  excess  of  the  pixel  dwell 
time  T. 

Physically,  Eq.  (12)  treats  the  clutter  return  as  being  unresolved  in  range 
coming  from  reflectors  located  L  m  in  front  of  the  radar.  The  £c(t)  process 
represents  mi cro-moti on ,  i.e.,  rapid  intra-pixel  variations,  of  the  clutter 
reflectors.  The  fc(t)  process  represents  macro-motion,  I.e.,  slow  gross-pixel 
movement,  of  the  clutter  reflectors.  The  combination  of  these  two  processes 

t 

£(t)  =  £<■  (t  )exp[j2*  /  fc(i)<fr].  (14) 

can  be  shown  to  have  a  spectrum 

S-~(  f )  ■  (2xo2)*1/Zexp[-(f-mfc)2/2a2],  (15) 

whose  .'ms  width  is 

o  ■  Uc2  +  °fc2)1/2-  U6) 

Thus,  the  long-term  bandwidth  of  our  clutter  model  is  a  composite  of  micro-motion 
and  macro-motion  c  ntributlons. 

3.0  PERFORMANCE  ANALYSES 

In  * h 1 s  section  we  shall  apply  the  models  presented  above  to  a  variety  of 
radar  -  .Tormance  calculations.  These  calculations  will  illustrate  the  interplay 
of  speckle,  frequency  instability,  and  clutter  effects. 
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3.1  VELOCITY  IHAoI.NO 

Consider  *  single-pixel  vel oc 1 ty- i meg i ng  radar  return  modeled  by  Eqs. 

( 1 ) - ( 4 ) ,  (s)-(lu),  (lib).  Let  us  assume  that  tne  target  range  L  is  known,  and 
that  a  surface  acoustic  wave  (SAW)  Fourier  transform  dev i ce/envel ope  detector 
combination8  is  used  to  obtain  the  processed  signal 

R(f)  -  |/dt  s*(t-2L/c)r(t)e*j2*ftJ 2.  (17) 

Typically,  R(f)  is  peak  detected,  i.e.,  f°  for  which  naxy  R(f)  • 

R ( f o )  s  ro  is  found;  the  radar  then  assigns  the  velocity  value 

•vz°  *  Xf°/2  and  the  intensity  value  R<>  to  this  pixel8.  Here,  we  shall 

study  the  velocity  resolution,  carrler-to-noi se  ratio  (CNR),  and  si gna I -to-nol se 

ratio  (SNR)  behavior10  of  the  processed  signal  R(f). 

The  mean  value  <R(f)>  of  the  processed  signal  consists  of  a  target-return 
contrl buti on 

<j/dt  £*(t-2L/c)£(t)e‘j2,ftj 2>  * 


EUV»)/d7|jCS,L)|  ‘•e'2aL25/Z  exp[-a(f-fQ)2/fRES2] 

*TfRES 


(18) 


where 


fRES  “  (25/Z/«T)[l+(*ofT)22L/ctf  +  T2/4tc2] 1/2,  (19) 

plus  a  noise  contribution 

<|/dt  s*(t-2L/c)n(t)e'jZ,ft|  z>  *  hvQ/n.  (20) 

Physically,  Tres  is  the  full  width  between  the  e"2  attenuation  points  of  the 
target-return  contribution  to  <R(f)>.  Tnus, 

VRES  “  XfRES/Z  (Z1) 


8.  Arsenault,  0.  R.,  and  Oolat,  V.  S. ,  "Compact  Multiple-Channel  SAW 

SI idinq-Window  Spectrum  Analyzer,"  1981  Ultrasonics  Symposium  (IEEE.  New 
York,  1981),  pp.  220-225. 


Biron,  0.  G.,  and  Edwards,  8.  E.,  "Moving  Target  Imaging  Radar  Utilizing 
Both  intensity  and  Velocity  Information,"  CLEO  ‘82  Technical  Digest  (Opt. 
Soc.  of  America,  Washington,  O.C.,  1982),  pp.  128-129. 

10.  As  elaborated  in  Ref.  1.,  CNR  measures  the  relative  strength  of  the  target 
return  and  noise  at  the  input  to  the  IF  signal  processor,  and  SNR  measures 
the  relative  strength  of  the  mean  target-return  and  fluctuations  at  the 
output  of  the  signal  processor. 
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is  the  velocity  resolution  of  the  radar  as  set  oy  tne  pixel  dwell  tine  (Tj, 
frequency  instaoility  ( 2 a f^L/ctf),  and  target  decorrelation  (tc) 
effects.  In  Fig.  3  we  nave  plotted  vpfg  vs. target  correlation  time  tc  for 
T  »  20  usee  (correspondi ng  to  dwell-time  limited  resolution  of  1.72  km/h  »  0.4s 
m/s),  assuming  no  frequency  instaoility  (2af2L/ctf  *  0).  This  figure  shows 
that  appreciable  resolution  loss  occurs  for  correlation  times  shorter  than  10 
nsec,  such  as  typically  occur  with  aerosol  targets.  In  Fig.  4  we  have  plotted 
vpEs  vs.  target  range  L  for  T*20  nsec,  and  four  pair  of  of,  tf  values, 
assuming  no  target  decorrel ati on  (tc  *  -) .  Here  we  see  that  the  dwell  time  Is 
sufficiently  short  that  severe  resolution  loss  occurs  only  under  extreme 
(<jf*200kHz,  tf  *  10" 4  sec)  instability  conditions. 

The  CNR  at  frequency  f  is  defined  to  be  the  ratio  of  the  signal  contribution 
to  <R(f)>  to  the  noise  contribution  to  <R(f)>.  At  the  target  Doppler  shift  f[) 
we  have  that 


CNR 


nE  (  \zp/ »)  Jdp[  p.L  ) 


v  e-2aL25/2 


"V^RES 


(22) 


Equation  (22)  can  be  recast  into  the  more  familiar  form6-11 


CNR 


nPT 


■Zal  Z 


5/2 


h  v  8 
o 


TIT 


RES 


(23) 


in  terms  of  the  transmitter  peak  power  ?j ,  an  effective  pulse  bandwidth 
B  =  Pj/E,  and  the  effective  heterodyne  detection  mixing  area  Ap.  From  ( 1 y ) 
we  see  the  final  fraction  on  the  right  in  (23)  represents  a  CNR  loss  that  is  due 
to  frequency  spreading  of  the  target  return.  This  loss  will  be  significant 
whenever  vrjs  is  appreciably  worse  than  the  1.72  km/h  dwell-time  limited 
value. 

The  SNR  at  frequency  f  is  defined  to  be 


SNR  •  (<R(f)>  -  hvQ/n)2/var(R(f)),  (24) 


i.e.,  it  is  the  ratio  of  the  squared  mean  target-return  contribution  to  R(f)  to 


11.  Shapiro,  J.  H.,  "Ta rget -Ref  1 ec t  i  v i ty  Theory  for  Coherent  Laser  Radars," 
App  1  .  Opt.  2_1 ,  3398-3407  (  1982  ). 
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;ne  varij'>.e  of  < . f  i ,  wnere  the  latter  Includes  target-return  and  noise 


fluctuations.  Equation  (24)  can  be  put  into  the  standard  form1 

CNR/2  _ 


SNR 


1  +  CNR/2SNR$at  +  (2CNR) 


-1 


in  terms  of  tne  CNK  and  the  saturation  si gna 1 -to-noi se  ratio 


SNRsat  =  <Ry(f)>2/var(Ry(f)), 


where 


(25) 


(26) 


R  (f )  •  j/dt  s*(t-2L/c)jt(t)  e*j2’ft|2  (27) 

is  the  target-return  component  of  R(f).  Equation  (25)  Implies  that  SNR  reaches 
the  limiting  value  SNRsat,  arising  from  target-return  fluctuations,  when 
CNR>>max( 5,2SNRsat ) .  After  some  tedious  algebra  we  find  that  the  saturation 
si gnal -to-noi se  ratio  at  frequency  fQ  obeys 

SNRsat  ‘  l2fR£S  /f,R£Sf"RES  *  ^  (28) 

where  f'ttES  1s  obtained  from  the  right  side  of  (19)  with  L«0  and  f "res  Is 
obtained  from  the  right  side  of  (19)  with  l  replaced  by  2L.  When  frequency 
instability  effects  are  insignificant  (as  Is  true  for  most  of  the  Fig.  4 
parameters),  (28)  reduces  to  the  well-known  unity  saturation  si  gnal -to-noi se 
ratio  of  laser  speckle.  However,  when  the  frequency  instability  term  dominates 
(19),  Eq.  (28)  yields  a  value  much  less  than  unity.  Indicating  worse  than  speckle 
target-return  fluctuations. 

3.2  RANGE  IMAGING 

Now  let  us  consider  a  single-pixel  range-imaging  radar  return  modeled  by 
Eqs.  (1),  (2),  (5),  (6),  (8)  -  (10),  (lib).  Here  we  assume  that  a  SAW  pulse 
compressor  is  used  to  obtain  the  processed  signal1* 

t(t)  »  j/dt  s*(t-2z/c)r(t) j 2.  (29) 

Typically,  i(z)  is  peak  detected,  l.e.,  z°  for  which  maxz  z(z)  *  t(z°)  = 
i°  is  found;  tne  radar  then  assigns  the  range  value  z°  and  the  Intensity 
value  to  this  pixel.  Once  again  we  shall  be  interested  in  the  resolution, 

CNR,  and  SNK  of  the  processed  signal;  the  development  parallels  that  just 
presented  for  velocity  Imaging. 
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The  average  target-return  contribution  to  the  processed  signal  is 


where 


<i  (j)»  *  /  dz‘  Q( 2  * )  F( z-z  * ;  2  ‘ ) . 
y  o 


4(2)  -  j / at  S* ( t-2z/ C )£( t ) | 2, 


Q ( 2 )  -  E( X2»1/Z/4) Jd7|j{p,z) | “cTe’*02  T(z) 


2  02 


is  the  infinite  resolution  average  target-returr  profile,  and 

2,1/2 


,2/,  2, 


F  (  2  ; 2  '  )  *  9(8/*2RES  )  exp  (-82  /2REs 


with 


cT/2 


1/2 


[1+(iWT/4) 2/(l+(»ofT) 222‘/ctf+T2/4tc2)]J 


(30) 


(31) 

(32) 


(33) 


(34) 


and 

g  -  /dzF ( z ; 2 '  )  * 

[l  +  (*ofT)  22z'/ctf+T2/4tc2+(trWT/4)  2]”1/Z,  (35) 

is  the  effective  range  probe  function  of  the  radar.  Equations  (30)  -  (34)  imply 
that  zRE 5  is  the  range  resolution  of  the  radar  at  range  z‘,  l.e.,  it  Is  the 
full  width  in  the  z  coordinate  between  the  e"2  attenuation  points  of  F(z;z'). 
Note  that  zp£5  includes  frequency  instability  and  target  decorrelation  effects 
In  addition  to  its  chirp  duration  and  chirp  bandwidth  H“pendence.  Under  typical 
operating  conditions,  the  target  range  profile  T(z)  may  be  a  preciably  narrower 
than  the  nominal  target  range  L.  In  this  case,  which  we  shall  employ  in  what 
follows,  we  may  set  z‘  *  l  1 n  (34)  and  (35).  To  exhibit  ti  ;  effect  of  target 
decorrelation  on  range  resolution,  we  have  plotted  zREs  vs*  lc  ’n  Fig.  5, 
with  T  *  10  usee  and  W  »  20  MHz  (corresponding  to  a  chirp-bandwidth  limited 
resolution  of  13.5m)  assuming  no  frequency  instability.  Here  we  see  the  same 
resolution  loss  that  transpired  for  the  Doppler  radar  in  Fig.  3,  viz.,  zrEs 
degrades  rapidly  for  tc  <10usec.  In  Fig.  6  we  have  plotted  zrEs  vs.  target 
range  L  for  T  *  10  usee,  U  ■  20  MHz,  and  three  pair  of  of,  tf  values, 


assuming  no  target  decorrelation,  ms  was  seen  in  F i y .  4  for  velocity  resolution, 
only  extreme  instaoility  conditions  cause  a  major  range  resolution  loss. 

The  CNR  at  range  z  is  defined  to  Be  the  ratio  of  <ty(z)>,  the  signal 
contribution  to  <i (z)>,  to  n»0/i,  the  noise  contribution  to  <i(z)>.  At  the 
nominal  target  range  l  we  have  that 


CNR  »  {  n / h  vQ ) /  d2*  g ( z  1  )  F( L-z '  ;L) 
0 


nPT  pAR  -2aL  . 

h!T5  .  2  nar  '* 

O  «L 


for  a  "hard"  target  obeying  (4), 


(nPT  3$  (  *  *  l )  z  '  res "R  -2  at.  , 

snri  ,Lz  son* 

for  a  "soft"  (aerosol)  target  obeying  (7), 

where  Py  Is  the  transmitter  peak  power,  B  =  Py/E  is  an  effective  pulse 
bandwidth,  Aps  Is  the  effective  heterodyne  mixing  area6,  and 

L  hard  *  9cT/21/2zRES, 


Z‘RES  »  cT»i/2/4[l+(.WT/4)z]1/2, 


^  soft  “  9cT’1/2/4z‘rES* 


Equations  (36a)  and  (37)  show  that  the  CNR  obtained  from  a  hard  target,  i.e.,  one 
that  is  unresolved  in  range,  takes  the  form  of  the  standard  monostatic  radar 
equation  times  the  k  factor  i-nard*  Via  (34)  and  (35)  it  is  easily  shown' 


lhard  *  Cl*(w*fT^2L/cVT‘/4tc2}'1/2’ 


which  Implies  fftard  with  ^nard  whenever  zpgs  is  approximately  equal 
to  its  chirp-band  width  limited  value.  The  hard  target  loss  factor  is  similar  to 
the  Doppler  radar  loss  factor  encountered  in  Eq.  (23);  somewhat  different 
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behavior  Is  seen  in  tne  soft  target  case.  Equation  (360)  is  Oastcally  toe 
usual  aerosol  monostatic  radar  equation  rimes  the  loss  factor  Lsoft,  »•*!  1  c h  we 
have  plotted  in  Fig.  7  vs.  tc  for  T  *  lUusec  and  W  *  20  MHz  assuming  no 
frequency  instaoility.  Comparing  Figs.  5  and  6,  we  see  that  for  aerosol  targets 
suostantlai  resolution  is  lost  before  any  perceptible  drop  in  CNR  occurs.  The 
reason  for  this  disparity  is  that  resolution  degrades  as  the  effective  probe 
function  F  broadens  due  to  loss  of  coherence.  Out  as  this  broadening  occurs  the 
radar  collects  aerosol  returns  from  a  larger  range  swath,  so  there  Is  almost  no 
immediate  reduction  in  CNR. 

The  SNR  at  range  z  is  defined  to  be  the  ratio  of  the  squared  mean 
target-return  contribution  to  £ ( 2 )  to  the  variance  of  i(z),  and  obeys  (25)  with 
CNR  from  (36)  and 

SNRsat  i  <ty(z)>2/var(  iy(z))  (41) 

giving  the  limiting  SNR  value  set  by  target-return  fluctuations.  We  shall  not 
give  explicit  results  for  SNRsat;  It  Is  sufficient  to  note  (cf.  (Z81)  tnat 
SNRsat  £  1  with  equality  when  frequency  instability  effects  are  Insignificant. 
3.3  CLUTTER-LIMITED  FALSE  ALARM  PROBABILITY 

In  this  section  we  shall  study  the  cl  utter- 1 i ml  ted  false  alarm  probability 
of  the  Doppler  radar  whose  single-pixel  vel oc 1 ty-i magi ng  performance  was 
addressed  in  Section  3.1.  Here  we  shall  consider  the  following  Idealized 
single-pixel  target  detection  problem1.  For  a  particular  range  L,  and  Doppler 
frequency  fg,  the  processed  signal  R ( f  q  )  defined  In  (17)  comprises  either 
clutter-return  plus  noise  (hypothesis  Hg)  or  target-return  plus  noise 
(hypothesis  Hi).  For  simplicity,  we  shall  assume  that  target  decorrelation  and 
frequency  Instability  effects  are  negligible  on  the  target  return,  and  frequency 
instability  effects  are  Insignificant  on  the  clutter  return.  The  radar  receiver 
will  use  thr  threshold  detection  rule 

say  H j 

«(f„)  7  Y  (42) 

say  Hq 

to  decide  on  the  absence  or  presence  of  a  moving  target  at  range  L  velocity 
Xfg/2  given  the  observation  R(fg).  Three  performance  quantities  of  interest 
are  the  no  1 se- 1 1  mi  ted  false  alarm  probability 
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•  Pr[K(fu)  >  f|£(t)  •  noise] 


(*3) 


tne  cl jtter-1 iml ted  false  alarm  probability 

PFc  *  PrC*( fQ)  1  f|i(0  *  Clutter  ♦  noise],  ( 44  ) 

ana  tne  aetectlon  probability 


p0  •  pfCR(fg)  2.  t (£.( t )  ■  target  ♦  noise].  (45) 

Ideally,  we  want  to  acnieve  Pp*  Ppc  «  0  and  Pq  ■  1,  but  because  of  noise 
0  <  Pp  and  P q  <  1  prevail,  ana  because  of  clutter  Pp  <  Ppc.  It  has  long 
been  known  for  the  problem  we  have  set  up  here  that1 

Pp  «  exp(-nT/hv0),  (46) 

and 


a  m  P  1/U*C*R) 


(47) 


where  CNR  Is  given  by  (23)  with  25/2/,Tfpps  •  i.  Thus,  for  example,  j  34  dB 
carrier- to-nolse  ratio  suffices  for  Pjj  «  0.99  operation  with  Pp  »  10- 12  In 
single-pixel  target  detection  In  the  absence  of  clutter.  When  clutter  Is  present 
and  HTI  Is  sought  at  velocities  near  zero  km/h  It  Is  Ppc  not  Pp  that  controls 
the  system's  false  alarm  rate. 

To  evaluate  Ppc  we  note  that  the  assumptions  we  have  made  Imply  that 
PfCR(^0)  >  t|.r(t)  ■  clutter  ♦  noise,  f  *  ^] 


1/Cl  ♦  CLNR(  f  rj  -  f  )] 


(48) 


where 


n£  /n  v 

CLMK(f)  *  - - - - -  exp 

[l*(«ocT)  Z/2]UZ 


(»fT/2)  2 
!♦(  »ocT)  2/ 2 


(*9) 


Is  the  c) utter-to-nol se  ratio  at  frequency  detuning  f.  Averaging  (48)  over  the 
fc  statistics  yields  tne  final  result 

PFc  •  /"«  <(2»0f2c)*1/2  exp[-(^a.fc)2/2of2c] 

l/[l*CLNR(f--*)] 

•PF  D  .  (50) 

U«  have  plotted  Ppc  vs.  CLNR,  the  cl utter-to-nol se  ratio  at  zero  detuning.  In 

Fig.  8  for  Pp  •  10'12,  T  ■  20  nsec ,  oc  -  30  kHz,  mpc  *  0  Hz, 

ofc  »  15  kHz,  and  two  values  of  fg  (corresponding  to  v2  »  vres/2  and  vz 

*  1.5  vre5,  where  vres  *  23/2X/*T  is  the  dwell-tlme  Halted  velocity 
resolution).  Fig.  8  shows  the  importance  of  clutter  when  trying  for 

h 1 gh-perforaance  HTI  at  velocities  comparable  to  the  ras  clutter  velocity  X(oc2 

♦  opc2)l/2/2. 

4.0  EXPERIMENTAL  RESULTS 

In  order  to  test  various  aspects  of  the  models  presented  In  Section  2.0,  we 
have  been  examining  the  statistics12  of  target  return  and  clutter  data  obtained 
with  the  MIT  Lincoln  Laboratory  Doppler  radar  test  bed9.  We  shall  report  herein 
some  of  our  preliminary  results. 


12.  Meslte,  P.  L.,  "Laser  Speckle  and  Clutter  Effects  on  Moving  Targets  Observed 
with  an  Optical  Radar,"  S.  M.  Thesis,  Dept,  of  Elect.  Eng.  ard  Comput. 

Scl.,  M.I.T.,  September  1983. 
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4.1  TAKUET  SPECKLE 

A  key  feature  of  tne  target-return  model  described  In  Section  2.1  is  tne 
presumption  that  tne  envelope  detected  return  from  a  moving  rough-surfaced  hard 
target  should  show  laser  speckle  fluctuations.  We  tested  this  presumption  by 
analyzing  returns  obtained  with  the  Doppler  test  bed,  from  a  truck-borne 
flame-sprayed  aluminum  calibration  plate  moving  at  constant  velocity.  Previous 
work3*11  with  the  M.I.T.  Lincoln  Laboratory  pul sed-lmager  test  bed13  had 
confirmed  the  diffuse  nature  of  this  calibration  plate  for  stationary  operation. 
Figure  9  shows  a  bar  chart  comprising  an  experimental  histogram  p(Y)  of  450 
spatially  sampled  normalized  speckle  plate  returns  Y  taken  by  the  Ooppler  radar 
with  the  target  at  2.7  km  range  moving  at  about  24  km/h.  Because  the  radar  uses 
a  linear  envelope  detector,  these  returns  should  be  Rayleigh  distributed  (with 
the  probability  density  shown  by  the  curve  In  Fig.  9)  If  our  assumed  speckle 
model  Is  correct3.  A  x2  goodness-of -f 1 t  test  for  the  Rayleigh  distribution  on 
this  histogram  gives  X2  *  9.426  with  11  degrees  of  freedom,  corresponding  to  a 
level  of  significance  between  0.5  and  0.75  and  excellent  agreement  between  theory 
and  data. 

The  Fig.  9  data  were  obtained  In  the  far  field  of  the  Doppler  radar. 

Similar  measurements  made  at  0.27  km,  which  Is  the  near  field  of  the  13  cm 
diameter  radar  telescope,  yielded  histograms  far  too  narrow  to  fit  a  Rayleigh 
distribution.  This  discrepancy  was  not  due  to  saturation  of  the  radar's 
electronics,  nor  could  It  be  explained  by  scanning-induced  decorrel atl on33 ;  it  Is 
a  subject  of  continuing  research. 

4.2  CLUTTER  SPECTRA 

The  essential  feature  of  the  clutter  model  given  In  Section  2.2  Is  Its 
distinction  between  micro-motion  and  macro-motion  contributions  to  the  frequency 
spectrum  of  clutter  returns.  To  test  this  decomposition,  we  analyzed33  clutter 
spectra  obtained  with  the  Doppler  test  bed  from  trees.  Figure  10  shows  an 
average  of  7920  Doppler  spectra  from  tree  clutter  at  0.5  km  range  obtained  under 
10-16  km/h  wind  conditions.  As  shown  In  this  figure,  the  central  peak  Is  well 
fit  by  the  Gaussian  profile  (15)  with  xifC  *  0  kHz  and  a  «  43.3  kHz;  the  region 
outside  of  the  central  peak  Is  ascribed  to  a  combination  of  SAW  sidelobes  and 
receiver  noise. 

13.  Harney,  R.  C. ,  and  Hull,  R.  J.,  “Compact  Infrared  Radar  Technology,"  Proc. 

SPIE  227,  162-170  (1980). 


F i 9 u r e  1U  represents  the  composite  of  micro-motion  and  macro-motion 
effects.  To  separate  these  components  In  the  data,  each  Individual  Doppler 
spectrum  was  windowed,  to  eliminate  outlying  points  due  to  sldelobes  and  noise, 
and  its  average  frequency  shift  computed.  The  windowed  spectra  were  then  shifted 
in  frequency,  so  that  their  average  frequency  shifts  all  fell  at  zero,  and 
averaged.  The  result  of  this  procedure,  shown  In  Fig.  11,  represents  the 
micro-motion  of  the  clutter.  As  can  be  seen  from  this  figure,  the  central  peak 
Is  well  fit  by  the  Gaussian  spectrum  (13)  with  oc  «  33.1  kHz.  Finally,  Fig.  12 
shows  the  histogram  of  average  frequency  shifts  for  the  7920  windowed  Doppler 
spectra.  This  histogram  represents  the  macro-motion  of  the  clutter,  and  Is 
compared  in  Fig.  12  with  the  Gaussian  density  (see  Section  2.2). 

PfcU)  *  (2*ofcar1/2  exp  [-U  -mfc  )  Z/2ofc  2]  (51) 

for  the  macro-motion  Instantaneous  frequency  shift.  He  see  from  Fig.  12  that 
(51)  with  mfc  *  34.9  kHz  and  afc  ■  35.6  kHz  gives  a  good  fit  to  the  data. 

Figures  11  and  12  are  not  fully  consistent  with  Fig.  10,  l.e.,  the  peak  In 
Fig.  10  should  be  fit  by  a  Gaussian  centered  at  mfC  «  34.9  kHz  with  rras 
bandwidth  (oc2  ♦  *  48.6  kHz  according  to  Figs.  11  and  12  and  the 

model  of  Section  2.2.  Although  the  above  rms  bandwidth  Is  not  too  far  from  the 
o  *  43.3  kHz  value  used  In  Fig,  10,  the  above  center  frequency  ca.inot  be  correct 
as  It  corresponds  to  27  m  of  net  tree  motion  away  from  the  radar  over  the  2.5 
minute  data  acquisition  period.  Given  the  33.7  kHz  spacing  between  adjacent 
frequency  samples  In  our  Doppler  spectra,  plus  the  problems  with  sldelobes, 
noise,  and  speckle  fluctuations,  however,  the  myc  discrepancy  Is  not  so 
unusua  1  . 

In  addition  to  examining  the  spectral  features  of  the  tree  clutter  data,  we 
have  also  considered12  the  amplitude  fluctuations  of  the  micro-motion 
contribution,  which  should  be  Rayleigh  distributed  according  to  the  model  In 
Section  2.2.  Preliminary  data  analysis  seems  to  Indicate  amplitude  fluctuations 
may  be  better  described  by  a  lognormal  rather  than  a  Rayleigh  distribution. 


5.0  CONCLUSIONS 
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Laser  radar  signal  models  were  presented  for  Doppler  Imagers  and 
ch 1 rped-pu 1 se  range  imagers.  These  models  were  used  to  assess  the  Impacts  of 
target  decorrelation  and  laser  frequency  Instability  on  the  velocity  and  range 
resolution  of  long  dwell  time  systems.  The  results  Indicated  that  decorrelation 
is  a  more  severe  proolem  than  frequency  Instability.  A  clutter  model,  containing 
micro-motion  and  macro-motion  components,  was  used  to  analyze  the  clutter-limited 
false  alarm  probability  In  single-pixel  detection.  The  appropriateness  of  the 
preceding  models  was  tested  using  data  from  the  M.I.T.  Lincoln  Laboratory  Doppler 
radar  test  bed.  Reasonable  confirmation  was  found  for  hard-target  speckle  In 
far-fleld  operation,  and  the  mlcro-motlon/macro-motlon  decomposition  of 
tree-clutter  spectra. 
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FIGURE  1.  SYSTEM  STRUCTURE  FOR  MONOSTATIC 
DOPPLER  OR  CHIRPED  LASER  RADAR. 


FIGURE  2.  NORMALIZED  SPECTRAL  DENSITIES 
FOR  THE  LASER  FREQUENCY- 
INSTABILITY  EFFECT  ON  THE  TAR¬ 
GET  RETURN:  (a)  the  very  short 
range  regime  2afz/c  <<  1,  given 
by  Eq.  (11a);  (b)  the  normal 
operating  regime  2z/ctf  <<  1, 
2a*z/c  >>  1,  given  by  Eq.  (lib); 
(c)  the  very  long  range  regime 
2z/ctf  >>  1,  given  by  Eq.  (11c). 


VELOCITY  RESOLUTION 

VMg  «  B^'X/irTXI  . 

X  ■  10. •  paa.  T  -  20  pa 


10"'  1  10  101 
1,  (pa) 


FIGURE  3.  VELOCITY  RESOLUTION  V-p-  vs. 

TARGET  CORRELATION  TIME  tc  FOR 
A  10.6  urn  WAVELENGTH  20  usee 
DWELL  TIME  DOPPLER  RADAR  WITH 
NO  FREQUENCY  INSTABILITY. 
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FIGURE  4.  VELOCITY  RESOLUTION  Vres  vs. 

RANGE  L  FOR  A  10.6  ym  WAVE¬ 
LENGTH,  20  ysec  DWELL  TIME 
DOPPLER  RADAR  VIEWING  A  NON- 
DKCORRELATING  TARGET. 
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RANGE  RESOLUTION  VS.  TAR-  FIGURE 

GET  CORRELATION  TIME  tc  FOR  A 
CHIRPED  LASER  RADAR  WITH  10  usee 
CHIRP  DURATION,  20  MHz  CHIRP 
BANDWIDTH,  AND  NO  FREQUENCY 
INSTABILITY. 


6.  RANGE  RESOLUTION  ZRES  vs .  RANGE 
L  FOR  A  CHIRPED  LASER  RADAR 
WITH  10  usee  CHIRP  DURATION  AND 
20  MHz  CHIRP  BANDWIDTH  VIEWING 
A  NON-DECORRELATING  TARGET. 


SOFT-TARGET  CNR  LOSS 

1— ,  «  (1  ♦  T*/4«//|1  ♦  (»WT/4rt  -*'*) 
T  •  10  M*.  *»  *  20  Mttl 


-»  -10  0  10  20  30  40  M 

CUM  |d» 


7.  SOFT  TARGET  CNR  LOSS  Lsoft  vs. 
TARGET  CORRELATION  TIME  tc  FOR 
A  CHIRPED  LASER  RADAR  WITH  10 
usee  CHIRP  DURATION,  20  MHz 
CHIRP  BANDWIDTH,  AND  NO  FRE¬ 
QUENCY  INSTABILITY. 


FIGURE  8.  CLUTTER-LIMITED  FALSE  ALARM 
PROBABILITY  PFc  vs.  CLUTTER- 
TO-NOISE  RATIO  AT  ZERO  DETUN¬ 
ING  CLNR  FOR  TWO  DOPPLER  SHIFT 
FREQUENCIES  fD;  Pp  =  10”12, 

T  =  20  usee,  ac  =  30  kHz,  mfc  ■ 
0  kHz,  and  ajc  *  15  kHz. 


CUmfA  INTtNUTY 


FIGURE  9.  HISTOGRAM  p(Y)  OF  TARGET  RETURN 
AilPLITUDES  i  FROM  A  MOVING 
SPECKLE  TARGET  (FLAME-SPRAYED 
ALUMINUM  CALIBRATION  PLATE) ; 
SMOOTH  CURVE  IS  THE  RAYLEIGH 
PROBABILITY  DENSITY. 


FIGURE  10.  AVERAGE  OF  7920  CLUTTER 

SPECTRA  (POINTS)  FROM  TREES 
AT  0.5  km  RANGE  IN  10-16 
km/ll  WIND  CONDITIONS  AND 
NORMALIZED  GAUSSIAN  (CURVE) 
FROM  (15)  WITH  m-  =  0  kHz 
AND  a  =  43.3  kHz  FIT  TO  CEN¬ 
TRAL  PEAK. 
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FIGURE  11.  AVERAGE  OF  7920  WINDOWED  CLUT¬ 
TER  MICRO-MOTION  SPECTRA 
(POINTS)  FROM  TREES  AT  0.5  km 
RANGE  IN  10-16  kra/h  WIND  CON¬ 
DITIONS,  AND  NORMALIZED  GAUS¬ 
SIAN  (CURVE)  FROM  (13)  WITH 
ac  -  33.1  kHz  FIT  TO  CENTRAL 
PEAK. 


FIGURE  12.  HISTOGRAM  p(fc)  OF  MACRO¬ 
MOTION  DOPPLER  SHIFTS  fc  OF 
7920  CLUTVER  SPECTRA  FROM 
TREES  AT  0.5  km  RANGE  IN  10- 
16  km/h  WIND  CONDITIONS: 
SMOOTH  CURVE  IS  GAUSSIAN 
PROBABILITY  DENSITY  WITH  mfc=> 
34.9  kHz  AND  afc  =*  35.6  kHz. 
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